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Foreword 

This  report  describes  results  obtained  in  the  second  year  of  a  three  year  research  study  on 
nonlinear  modeling  and  control  of  flexible  space  structures  with  application  to  rapid  slewing 
and  precision  pointing  of  space-based,  directed  energy  weapons.  The  project  is  funded  by 
SDIO/IST  and  managed  by  AFOSR/SDIO  (AFSC).  Results  reported  herein  are  for  the 
period  1  Sept.  1988  -  31  Aug.  1989. 

The  project  is  managed  by  Ltn.  Colonel  .J.  Crowley  and  Dr.  A.  Amos.  We  wish  to  thank 
both  of  these  individuals  for  their  insight  and  direction  on  this  project. 


Aeeesslon  For 

NTIS  CRAil 

DTIC  TAB 

□ 

Unannounced 

□ 

Just  ii’icntion _ 

By 

Distribution/ 


Availability  Coden 
Avail  and/or 
Special 


Dlst 


4 


CONTENTS  ii 

Contents 

1  Research  Objectives  and  Project  Summary  1 

2  Feedback  Linearization  and  Stabilization  of  Nonlinear  Systems  5 

2.1  Computation  of  Partial  Linearizing  Feedback  Compensation .  7 

2.2  Nonlinear  System  Transmission  Zeros .  9 

2.3  PLF  Computations  for  Nonlinear  MIMO  Systems .  10 

2.4  Partial  Linearization  and  Variable  Structure  Control  Systems .  13 

2.5  Robust  Stabilization  of  Nonlinear  Systems .  16 

2.6  Partial  Linearization  for  Lagrangian  Systems .  19 

3  Reduced  Order  Modeling  of  Multibody  Systems  and  Nonlinear  Control 

Design  22 

3.1  PFL  with  Reduced  Order  Models .  24 

3.1.1  Quasi- Rigid  Model .  26 

3.1.2  The  reduced  flexible  model .  28 

3.1.3  Feedback  linearization  and  decoupling  of  the  quasi-rigid  model  ....  29 

3.1.4  Partial  feedback  linearization  of  the  Reduced  Flexible  Model .  32 

3.2  Conclusions  and  Directions  .  32 

4  Some  Design  Approaches  for  Combined  Rapid  Slewing  and  Precision  Point¬ 
ing  34 

4.1  Practical  Implementation  of  Time-optimal  Slewing  for  an  Inertial  Load  ...  37 

4.2  Implementation  of  Rapid  Slewing  with  Coordinated  Modes  of  Actuation  .  .  41 

4.3  Coordinating  Control  for  System  Slewing  and  Multibody  Alignment .  46 

5  Simulation  Tradeoff  Studies  for  Multiaxis  Slewing  of  SBL  System  Model  49 

5.1  Multibody  System  Model  for  Multiaxis  Slewing  and  Precision  Alignment  .  .  49 

5.1.1  FEM:  Collocation  by  Splines .  51 

5.1.2  Reduction  of  the  Kinetic  Energy  Function .  51 

5.1.3  Reduction  of  the  Potential  Energy  and  Dissipation  Functions .  52 

5.1.4  Lagrange’s  Equations .  54 

5.2  Computer  Simulation  Model  for  Nonlinear  Control  Law  Tradeoff  Studies  .  .  54 

5.2.1  Simulation  model  geometry  assumptions  and  parameters .  55 

5.3  Simulation  Results  for  PLF  Control  with  Rapid  Slewing  and  Precision  Pointing  56 

6  Conclusions  and  Directions  61 

A  Advantages  of  Gibbs  vector  description  of  nonlinear  kinematics.  65 

B  Supporting  Computations  for  area  moment  of  inertia  and  B-spline  model 
for  3-axis  slewing  with  tripod  shaped  appendage  67 

C  Simulated  Time  Responses  for  Multiaxis  Slewing  73 


LIST  OF  FIGURES  iii 

List  of  Figures 

2.1  Nonlinear  Control  Concept  Using  Dynamical  System  Inverse .  7 

2.2  Partial  Linearizing  Feedback  via  Nonlinear  Inverse  Model .  11 

2.3  Partial  Feedback  Linearization  and  Zero  Dynamics  .  11 

4.1  Implemenation  of  Time-Optimal  Servo .  41 

4.2  Root  Locus  for  Linear  Mode  Sensitivity  of  Saturation  Mode  Slewing  ....  41 

4.3  Coordinated  Rapid  Slewing  Control  Implementation .  46 

5.1  Simulation  model .  55 

2.1  Cross-section  of  tripod  metering  truss .  67 


List  of  Tables 


5.1  Standard  Notation  for  Lagrangian  Mechanics .  50 

5.2  Matrix  of  Simulation  Models  Considered .  57 

5.3  Matrix  of  PLF/Slewing  Control  Implementations  .  57 

5.4  Slewing  Times  for  Standard  3-axis  Maneuver .  58 

5.5  Relative  Peak  Torque  Increase  with  PLF  Compensation  DOF  .  59 


TSI-89-12-12-WB 


1 


1  Research  Objectives  and  Project  Summary 

The  objectives  of  the  research  project  are  to  develop  nonlinear  control  design  techniques 
based  on  the  idea  of  feedback  linearization  for  application  to  control  of  flexible  space  struc¬ 
tures.  The  first  year  effort  focused  on  the  development  of  a  generic  model  of  the  dynamics  of 
a  multi-body  system  with  elastic  interactions  undergoing  large  angle  slewing  motions.  The 
model  was  specialized  and  scaled  to  represent  available  data  on  the  SBL  prototype  system 
and  a  computer  simulation  was  developed.  The  significance  of  partial  feedback  linearizing 
compensation  where  principal  nonlinear  couplings  are  dynamically  compensated  for  by  non¬ 
linear  feedback  control  was  demonstrated  for  the  problem  of  rapid  and  precision  large  angle 
slewing  of  a  central  rigid  body  with  elastic  structure.  To  achieve  large  angle  slewing  de¬ 
coupling  control  was  used  in  concert  with  nonlinear,  time-optimal  switching  control  for  the 
rapid  slew.  Results  were  presented  at  the  27th  IEEE  C'DC  in  Austin,  Texas. 

At  the  request  of  Dr.  A.  Amos,  second  year  project,  activities  included  the  development, 
of  a  detailed  simulation  model  including  multiaxis  attitude  dynamics  and  structural  flexure. 
The  model  is  described  in  detail  in  Section  5.1  along  with  parameters  chosen  for  simulation 
tradeoff  studies.  The  system  model  parameters  were  chosen  from  available  data  to  approx¬ 
imately  model  the  dynamics  of  the  test  article  under  construction  at  the  ASTREX  facility 
at  Air  Force  Astronautics  Laboratory. 

Implementation  Alternatives  for  Feedback  Linearizing  Control.  In  the  second  year 
we  have  focused  on  several  practical  aspects  of  the  implementation  of  such  decoupling  control 
laws.  Reduced  order  modeling  of  the  flexible  structure  elastic  dynamics  was  considered  for 
implementation  of  nonlinear  feedback  compensation.  A  time  scale  analysis  was  performed 
and  simulation  studies  were  performed  to  illustrate  tradeoffs  in  LOS  slewing  and  pointing 
precision  vs.  peak  torque  requirements.  A  relatively  soft  structural  stiffness  model  was  used 
to  illustrate  the  tradeoffs.  The  time  scale  analysis  based  on  singulc.r  perturbations  identified 
several  alternatives  for  enhanced  precision  of  decoupling  control  implementation  without 
model  order  increase.  The  use  of  structural  actuators  for  deformation  shaping  in  concert 
with  slewing  control  was  considered  and  additional  simulation  studies  are  being  performed. 
Considerations  for  robust  implementation  of  slewing  control  has  focused  on  precision  of 
the  decoupling  compensation  and  stability  robustness  with  parasitic  dynamics  arising  from 
elastic  structure  response.  This  question  directs  attention  to  the  sensitivity  of  the  control 
to  structural  damping  and  will  provide  direction  for  experiments  to  be  designed  in  the  last 
quarter  of  this  years  effort. 

A  critical  observation  in  our  studies  of  control  architectures  for  SBL  systems  is  the  inte¬ 
gration  of  a  variety  of  actuators  for  spacecraft  attitude  control  (e.g.  thruster  jets,  moment  um 
wheels,  CMG’s,  etc.),  multibody  articulation,  optical  system  components  (e.g.  steering  and 
deformable  focusing  mirrors),  and  structural  vibration  control  (e.g.  proof  mass  devices, 
embedded  piezoelectrics,  etc.)  to  achieve  principal  system  performance  objectives.  In  this 
report  we  describe  alternative  implementations  for  coordinating  control  activity  between 
actuators  of  different  types. 
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Tradeoff  Studies  in  Implementation  of  Nonlinear  Decoupling  Control  Performed. 

The  idea  behind  feedback  linearizing  compensation  is  to  compensate  for  nonlinear  couplings 
by  cancelation  of  critical  nonlinear  terms  effecting  system  dynamics.  Implementation  in¬ 
volves  the  introduction  of  additional  control  authority  and  may  suggest  the  use  of  special 
actuator  configurations.  Coordinating  control  activity  between  various  types  of  actuators  is 
an  important  feature  of  such  control  laws.  As  part  of  the  second  year  effort  detailed  simula¬ 
tion  studies  were  performed  to  delineate  tradeoffs  in  implementation  of  feedback  linearizing 
control  of  a  SBL  system  model  undergoing  multiaxis  slewing. 

An  important,  alternative  implementation  is  possible  based  on  the  use  of  fast,  switching 
control  laws.  The  implementation  of  feedback  linearization  by  smooth  control  involves  ex¬ 
plicit  compensation  for  critical  nonlinear  dynamic  couplings.  In  the  current  study  we  have 
shown  that  feedback  linearization  can  be  achieved  by  implicit  compensation  based  on  the 
notion  of  variable  structure  or  sliding  mode  control  (see  Section  2.4).  In  this  case  the  control 
law  is  discontinuous  with  respect  to  a  switching  or  sliding  surface  in  the  model  (slow  time 
scale)  state  space.  An  important  advantage  is  that  the  control  law  can  be  implemented  based 
only  on  measurements  of  primary  body  attitude  parameters.  Given  sufficient  control  author¬ 
ity  to  maintain  sliding  robust  performance  can  be  maintained  independent  of  variations  in 
flexible  multibody  dynamics. 

New  Results  Obtained.  A  potentially  important  feature  of  feedback  linearization  is  the 
extent  to  which  the  idea  can  be  integrated  with  other  standard  design  methods  for  multi- 
variable  control  systems.  We  have  described  several  alternatives  which  indicate  advantages 
of  the  approach  for  control  of  flexible  structures.  In  particular,  the  integration  of  relatively 
weak  structural  actuation  for  vibration  suppression  can  be  incorporated  with  the  approach. 
The  idea  of  deformation  shaping  control  has  been  developed  by  Dr.  T.A.W.  Dwyer  is  re¬ 
ported  in  several  papers  prepared  as  part  of  this  project.  An  alternative  approach  which  is 
briefly  described  in  this  report  is  to  study  the  tradeoff  between  implementation  costs  for  feed¬ 
back  linearizing  control  with  passive  structural  damping  vs.  active  structural  damping.  An 
approach  to  the  required  feedforward  compensation  is  developed  based  on  computations  of 
Partial  Linearizing  Feedback  compensation  in  terms  of  the  normal  form  equations  of  Byrnes 
and  Isidori. 

Robust  s'abilization  and  control  performance  of  nonlinear  systems  is  an  important  issue 
in  applications.  Several  techniques  for  robust  stabilization  have  been  developed  by  Spong 
and  Vidyassagor,  Slotine  and  Sastry,  C'orless,  Gutman,  and  others.  These  methods  all  rely  on 
special  structure  matching  conditions  for  the  model  uncertainty  in  developing  robust  control 
design  methods.  In  this  research  we  have  developed  a  new  approach  which  utilizes  the  frame¬ 
work  of  adaptive  control  to  provide  robust  stabilization  of  nonlinear  systems.  The  approach 
does  not  rely  on  structure  matching  conditions  and  represents  a  significant  improvement  over 
available  methods. 

Professional  Personnel  The  principal  investigator  for  this  project  is  Dr.  William  H. 
Bennett  and  co-principal  investigators  are  Drs.  Harry  G.  Kwatny  and  Gilmer  L.  Blankenship 
from  TSI.  Dr.  Thomas  A.  W.  Dwyer  was  consultant  on  the  project.  We  would  also  like  to 
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acknowledge  the  parttime  support  from  Dr.  Oussima  Akhrif  who  graduated  in  July  1989 
with  Ph.  D.  from  Electrical  Engineering  Department  of  the  University  of  Maryland.  Dr. 
Akhrif's  dissertation  developed  aspects  of  feedback  linearization  for  the  multibody  SBL 
models  developed  in  this  study  and  some  important  results  are  summarized  in  this  report. 

Technical  Reports/Presentations.  As  part  of  the  research  program  we  have  organized 
two  invited  technical  sessions,  presented  several  technical  papers,  and  submitted  several 
additional  papers  for  publication  as  follows. 

1.  Nonlinear  Dynamics  and  Control  of  Aerospace  Systems ,  invited  session  at  27t.h  IEEE 
Cntrl.  Dec.  Conf.,  Austin,  TX,  Dec.  1989. 

2.  Robust  Control  of  Uncertain  Nonlinear  Systems,  invited  session  at  1989  Amer.  Cntrl. 
Conf.,  Pittsburgh,  PA,  June  1989. 

Publications/Presentations 

1.  H.G.  Kwatny  and  W.H.  Bennett,  “Nonlinear  Dynamics  and  Control  Issues  for  Flexible 
Space  Platforms,”  Proc.  IEEE  Cntrl.  Dec.  Conf.,  Austin,  TX,  Dec.  1988. 

2.  T.A.W.  Dwyer,  III,  “Slew-Induced  Deformation  Shaping”,  Proc.  IEEE  Cntrl.  Dec. 
Conf..  Austin,  TX,  Dec.  1988. 

3.  0.  Akhrif,  G.  L.  Blankenship,  and  W.H.  Bennett,  “Robust  Control  for  Rapid  Reorien¬ 
tation  of  Flexible  Structures,”  Proc.  1989  Amer.  Cntrl.  Conf.,  Pittsburgh,  PA,  June 
1989. 

4.  H.G.  Kwatny  and  H.  Kim,  “Variable  Structure  Control  of  Partially  Linearizable  Dy¬ 
namics,”  Proc.  1989  Amer.  Cntrl.  Conf.,  Pittsburgh,  PA,  June  1989. 

5.  W.H.  Bennett,  “Frequency  Response  Modeling  and  Control  of  Flexible  Structures: 
Computational  Methods,”  3rd  Annual  Conf.  on  Aerospace  Computational  Control, 
Oxnard,  CA,  Aug. 

6.  T.A.W.  Dwyer,  III  and  F.K.  Kim,  “Nonlinear  robust  Variable  Structure  Control  of 
Pointing  and  Tracking  with  Operator  Spline  Estimation”,  Proc.  IEEE  International 
Symposium  on  Circuits  and  Systems,  Covaliis,  Oregon,  May  9-11,  1989,  Paper  No. 
SSP15-5. 

7.  T.A.W.  Dwyer,  III  and  F.K.  Kim,  “Bilinear  Modeling  and  Estimation  of  Slew-Induced 
Deformations,”  J.  Astro.  Sci.  submitted. 

8.  T.A.W.  Dwyer,  III,  “Slew-Induced  Deformation  Shaping  on  Slow  Integral  Manifolds,” 
Control  Theory  and  Multibody  Dynamics,  Eds.  J.  Marsden  and  P  S.  Krishnaprasad, 
(to  appear),  Amer.  Math.  Society. 
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9.  T.A.W.  Dwyer,  III,  F.  Karray,  and  W.H.  Bennett,  “Bilinear  Modeling  and  Nonlinear 
Estimation”,  Proc.  Flight  Mechanics /Estimation  Theory  Symposium,  NASA  Goddard 
Space  Flight  Center,  May  1989. 

10.  T.A.W.  Dwyer,  III  and  J.  R.  Hoyle,  Jr.,  “Elastically  Coupled  Precision  Pointing  by 
Slew-Induced  Deformation  Shaping,”  Proc.  1989  Amer.  Cntrl.  Conf.,  Pittsburgh,  PA, 
June  21-23,  1989. 

11.  T.A.W.  Dwyer,  III  and  Jinho  Kim,  “Bandwidth-Limited  Robust  Nonlinear  Sliding 
Control  of  Pointing  and  Tracking  Maneuvers,”  Proc.  1989  Amer.  Cntrl.  Conf.,  Pitts¬ 
burgh,  PA,  June  21-23,  1989. 

12.  T.A.W.  Dwyer,  III  ,  F.  Karray  and  Jinho  Kim,  “Sliding  Control  of  Pointing  and  Track¬ 
ing  with  Operator  Spline  Estimation,”  3rd  Annual  Conf.  on  Aerospace  Computational 
Control,  Oxnard,  CA,  Aug.  28-30,  1989. 
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2  Feedback  Linearization  and  Stabilization  of  Nonlinear  Systems 

Conventional  techniques  for  stabilization  of  nonlinear  systems  via  feedback  control  are  still 
very  limited  and  tend  to  be  tailored  to  specific  situations.  Among  the  most  promising  new, 
general  approaches  utilize  linearization  (local  or  possibly  global)  by  Exact  Feedback  Lin¬ 
earization  (EFL)  [HSM83,  KC87).  EFL  methods  are  based  on  earlier  work  of  Krener  [Kre73] 
and  Brockett  [Bro78]  which  demonstrated  that  a  large  class  of  nonlinear  dynamical  systems 
can  be  exactly  (i.e.  globally  linearized)  by  a  combination  of  nonlinear  transformation  of 
the  state  coordinates  with  nonlinear  state  feedback.  More  recently,  the  connection  between 
these  methods  and  the  idea  of  input-output,  (or  Partial)  Linearizing  Feedback  (PLF)  by  con¬ 
struction  of  a  system  inverse  [Hir79]  has  been  articulated  in  a  series  of  papers  by  Byrnes 
and  Isidori  [BI85,  BI84] .  These  connections  have  engendered  a  series  of  design  methods  with 
representative  results  for  specific  applications  by  Kravaris  and  Chung  [KC87]  and  Fernan¬ 
dez  and  Hedrick  [FH87].  In  this  section  we  will  show  how  fundamental  these  constructions 
can  become  in  control  system  design,  discuss  alternatives  for  implementation,  and  suggest 
some  approaches  to  integrating  the  nonlinear  design  philosophy  with  more  conventional  ap¬ 
proaches.  We  focus  attention  in  this  section  on  fundamental  concepts  culminating  in  the 
description  of  the  design  approach  for  multibody  systems  from  the  perspective  of  Lagranian 
mechanics. 

The  idea  behind  feedback  linearization  is  conceptually  simple.  We  start  with  a  nonlinear 
system  model, 

*  =  f(x)  +  G{x)u,  (2.1) 

y  =  h(x),  (2.2) 

where  x  £  5Rn,  u,  y  £  3?m  with  G  —  [<71 , .  . . ,  gm]  and  assuming  the  vector  fields  /,  gi  are  (7°° 
for  each  i  =  1, .  . .  ,m  and  /( 0)  =  0.  The  model  structure  assumes  that  the  control  u  enters 
linearly.  The  feedback  linearization  problem  is  to  find  a  change  of  basis  in  the  state  space, 
z  —  T(x ),  with  T  diffeomorphic  and  a  feedback  law, 

u  =  a(x)  +  [3(  x)v, 

such  that  in  the  new  (z,v)  coordinates  the  (closed  loop)  model  has  the  form, 

z  =  Az  +  Bv. 

We  remark  that  if  it  possible  to  find  such  a  control  law  then  the  linearization  is  achieved 
through  the  introduction  of  active  control  authority.  An  important  feature  for  control  system 
design  is  that  the  range  of  validity  of  the  linearization  is  given  by  the  transformations, 
T(x),a(x),f3(x).  The  functions  may  be  defined  locally  or  globally. 

In  contrast,  the  conventional  approach  to  control  design  would  be  based  on  a  linear  model 
obtained  by  Taylor  expansion  of  the  vector  fields  about  given  equilibrium  conditions ;  xeq,ueq, 
satisfying, 

0  =  /(**,)  +  G(xeq)ueq. 
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The  conventional  linear  model  used  for  control  design  represents  perturbation  dynamics  with 
respect  to  the  equilibrium  conditions  and  assumes  the  form, 

(Ar)  =  A  Ax  -f  B  A u, 


where  Ax 


y  =  C  Ax 


x  —  xeq ,  Au  =  u  —  uea  and 


<•<? 


■  B  =  C(.r)l  , 


A  major  source  of  model  uncertainty  in  linear  control  system  design  arises  from  assumptions 
leading  to  the  linear  perturbation  model.  In  many  cases  it  may  be  difficult  to  estimate  the 
domain  of  attraction  for  the  equilibria.  Indeed,  in  aerospace  applications  the  control  design  is 
often  based  on  a  combination  of  gain  scheduling  to  take  into  account  the  dependence  of  linear 
perturbation  models  on  operating  point  conditions  which  are  subject  to  variation  (e.g.  trim 
conditions  in  aircraft  flight  control.)  For  example,  the  function  of  a  conventional  autopilot 
for  aircraft  is  to  compensate  for  changes  in  trim  conditions  and  provide  stabilization  so  that 
the  pilot  “feels'1  a  standard,  linear  response  to  stick  commands. 

Although  the  concept  of  feedback  linearization  in  control  system  design  is  potentially 
revolutionary,  its  application  has  many  antecedents  in  applications.  The  significance  for 
nonlinear  control  of  flexible  space  structures  is  the  emerging  technology  for  active  control 
and  sensing,  the  dynamics  associated  with  the  CSI  technology,  and  the  ability  of  a  com¬ 
prehensive  approach  to  nonlinear  dynamic  modeling  and  control  design  offered  by  the  ap¬ 
proaches  discussed  in  this  report.  Feedback  linearization  functions  in  certain  applications 
in  a  manner  similar  to  gain  scheduling  [MC80],  however,  linearization  is  achieved  about  a 
“nominal  model11  rather  than  about  an  operating  point.  Thus  equilibria  conditions  do  not 
arise  explicitly  in  linearization.  One  view  of  such  a  controller  structure  is  illustrated  in  block 
diagram  of  Figure  2.1.  The  process  linearization  which  facilitates  the  design  of  the  linear 
controller  is  obtained  by  the  introduction  of  an  Inverse  Force  Model  (IFM)  for  the  nonlinear 
multibody  system.  The  inverse  force  model  transforms  commanded  accelerations,  ac,  into 
equivalent  system  generalized  forces,  /.  Thus  the  linear  controller  is  designed  to  yield  desired 
system  accelerations  given  the  generalized  coordinates,  q,  and  their  rates,  q. 

Precursors  to  the  idea  of  feedback  linearization  is  pervasive  in  control  applications  abound. 
With  the  development  of  the  geometric  theory  of  nonlinear  systems,  computational  tools  and 
design  methods  are  becoming  available  to  address  control  system  design  on  a  much  larger 
scope.  It  is  clear  that  the  concept  of  feedback  linearization  is  pervasive  in  many  fundamental 
control  methods.  Our  study  has  focused  on  the  considerations  for  practical  implementation 
of  feedback  linearization  for  rapid  slewing  and  precision  pointing  of  aerospace  systems  with 
multibody  and  elastic  interactions.  Implementation  of  feedback  linearization  definitely  re¬ 
quires  enhanced  control  authority  and  issues  related  to  technology  for  control  actuation  will 
either  enable  or  restrict  its  application. 

It  is  often  suggested  that  controllers  designed  based  on  feedback  linearization  may  be 
sensitive  to  model  assumptions  since  the  linearization  is  achieved  by  cancelation  of  cer¬ 
tain  nonlinear  terms  in  the  system  model.  Our  study  of  rapid  slewing  control  of  a  generic 
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Figure  2.1:  Nonlinear  Control  Concept  Using  Dynamical  System  Inverse 

space-based  laser  model  has  shown  that  modeling  sensitivity  and  robustness  can  be  obtained 
through  judicious  application  of  control  authority.  A  central  issue  in  nonlinear  control  design 
is  the  limits  of  available  control  authority.  For  example,  actuator  saturation  contributes  to 
limits  on  control  authority.  Discontinuous  or  saturation  mode  operation  of  control  actua¬ 
tors  can  often  be  desirable  but  such  considerations  are  often  not  addressed  by  linear  design 
methods.  We  have  shown  that  specific  consideration  for  coordination  of  discontinuous  and 
continuous  modes  of  control  actuation  in  slewing  control  methods  for  multibody  systems 
can  be  readily  found.  Finally,  it  is  clear  that  the  cost  of  feedback  linearization,  in  terms 
of  increased  control  authority,  sensor  measurement  complexity,  or  computational  burden  for 
online  implementation  may  not  always  be  necessary  to  achieve  system  performance  objec¬ 
tives.  We  have  also  demonstrated  that  such  methods  can  be  readily  integrated  with  standard 
approaches  for  linear  design  once  a  primary  system  control  objective  is  identified  in  terms 
of  a  primary  system  output.  As  we  shall  show  the  role  of  conventional  linearization  and 
control  design  can  be  relegated  to  a  subsystem  whose  dynamics  are  decoupled  (by  the  action 
of  PLF)  from  a  set  of  system  primary  outputs. 

2.1  Computation  of  Partial  Linearizing  Feedback  Compensation 

Partial  linearization  derives  directly  from  the  Byrnes-Isidori  normal  form  for  nonlinear  sys¬ 
tems.  The  essentials  of  the  approach  are  most  easily  developed  for  single-input,  single  output 
systems  and  we  will  present  the  approach  in  that  context.  The  theory  for  extending  these 
results  for  multi-input,  multi-output  problems  is  now  complete  and  references  are  included. 
Consider  a  nonlinear  dynamical  system  in  the  form, 

x  =  f{x)  +  g(x)u  (2.3) 

y  =  h{x)  (2.4) 

where  /,  g  are  smooth  C'°°  vecter  fields  on  ??n  and  h  is  a  smooth  function  mapping  9?n  — > ■  9?. 
Now  if  we  differentiate  (2.4)  wt  obtain 

dh 

V  =  £-( /(* )  +  9(x)u)- 
ox 


(2.5) 
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In  the  case  that  the  scalar  coefficient  of  u  (viz.  |^(:r))  is  zero  we  can  differentiate  again  until 
a  nonzero  control  coefficient  appears.  The  number  of  required  differentiations  is  fundamental 
system  invariant  which  plays  a  role  in  constructing  a  system  inverse  and  therefore  in  PLF. 
The  Byrnes-Isidori  analysis  shows  that  this  integer  number  is  analogous  to  the  relative  degree 
for  a  linear  system  [BI84] . 

The  above  construction  can  be  made  precise  using  the  notation  of  differential  geometry 
which  has  found  application  in  analytical  mechanics  [Arn78j.  We  will  need  only  the  notion 
of  Lie  derivative  and  Lie  bracket.  The  Lie  (directional)  derivative  of  the  scalar  function  h 
with  respect  to  the  vector  field  /  is 

I,(/>)  =  (2.6) 

Since  the  above  operation  results  in  a  scalar  function  on  5?n,  higher  order  derivatives  can  be 
successively  defined 

L)(h)=Li(l)-'{,h))-.=  (il)-'(h),j).  (2.7) 

Then  we  can  write  (2.5)  as 

V  =  ( dh,f)  +  (dh,g)u 

=  Lf(h)  4-  Lg{h)u.  (2.8) 

If  Lg(h )  =  0  then  we  differentiate  again  to  obtain 

y  =  (dLf(h),  f)  +  [dlf(h),g)u 

=  L2f(h)  + Lg(Lf(h))u.  (2.9) 

If  Lg(Lkf~1(h ))  =  0  for  k  =  1, . . .  ,r  —  1,  but  Lg{LTf~l{h))  ^  0  then  the  process  terminates 
with 

^=L',(h)  +  L,[L'f'(h))u.  (2.10) 

The  system  (2.10)  can  be  effectively  inverted  by  introducing  a  feedback  transformation  of 
the  form 

u  =  i,(irWv  ~ LW]  <2'n) 

which  results  in  an  input-output  response  from  v  — ►  y  given  by 


a  linear  system. 

The  integer  r  >  0  can  be  viewed  as  a  relative  degree  for  the  nonlinear  system  (2.3)-(2.4). 
Note  that  if  we  define  new  state  coordinates  :  £  3?r  as 


Zk  =  L)-'(h),  k=  1 . i 


(2.12) 
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for  the  r-dimensional  nonlinear  system  (2.10),  then  the  system  model  can  be  written  in  state 
space  form  as, 

0 

0  1  0  •••  0  ' 


0  0  1 
0  0  0 
0  0  0 


0 

1 

0 


z  + 


0 


a(x)  +  p(x)u  J 


(2.13) 


where 


a(x)  =  L}(h),  p(x)  =  Lg{Uf\h)). 


(2.14) 


More  generally,  using  the  new  coordinates  c  (2.12)  and  introducing  a  nonlinear  feedback 
control  of  the  form 

(v  —  <t(x)) 

v  v  ”  (2.15) 


where 


u  = 


r  —  1 


p{x) 


a(x)  =  Y,faL)(h)  +  Lrf(h), 


k= 0 


p(x)  =  Lg(Lrfl(h)), 


(2.16) 

(2.17) 


with  0k  for  k  =  0, . . . ,  r  —  1  real  positive  coefficients  then  the  equations  (2.8)-(2.9)  can  be 
written  in  ‘reduced’  form; 


■  0 

1 

0 

0  ■ 

( 0  \ 

0 

0 

1 

0 

0 

z  — 

0 

0 

0  ••• 

0 

3  + 

.  -0o 

-01 

-07  ... 

0r  —  \  . 

1 1 ) 

y  =  [1,0 —  ,0]  r. 


(218) 

(2.19) 


2.2  Nonlinear  System  Transmission  Zeros 

Note  that  the  process  leading  to  (2. 18)— (2. 19)  provides  an  equivalent  state  space  realization 
for  the  v  >-+  y  input-output  response  of  McMillan  degree  r  <  n  (the  dimension  of  the  original 
state  space  model  (2.3)— (2.4))  by  decoupling  a  portion  of  the  system  dynamics  from  the  output 
response.  This  is  depicted  in  Figure  2.3.  Thus  the  new  state  coordinates  z  are  a  ‘partial’  state 
for  the  system.  Thus  stabilization  of  ( 2. 18)— (2. 19)  cannot  guarantee  stabilization  of  the  full 
state  model  (2.3)-(2.4).  We  remark  that  in  the  case  that  h(x)  is  such  that  the  relative  degree 
r  =  n  then  the  state  space  transformation  (2.12)  for  k  =  1, . . . ,  r  together  with  the  feedback 
transformation  (2.15)  exactly  linearizes  the  full  system  state  model  (2.3).  The  methods 
described  in  [HSM83]  identify  necessary  and  sufficient  conditions  for  the  existence  of  a  C°° 
function  h(x)  such  that  r  =  n  and  provides  a  computational  approach.  The  necessary  and 
sufficient  conditions  for  (local)  EFL  are  nongeneric  and  not  likely  to  be  satisfied  in  general. 
In  the  sequel  we  show  how  that  essential  property  of  involutivity  of  the  /  and  g  vector 
fields  will  almost  never  be  satisfied  for  realistic  models  of  flexible  space  structures  due  to  the 
infinite  dimensional  nature  of  the  state  space. 
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Byrnes  and  Isidori  [BI85j  describe  the  transformation  of  (2.3)-(2.4)  to  a  normal  form  in 
which  the  feasibility  of  PLF  control  can  be  assessed.  The  main  result  provides  the  existence 
of  a  diffeomorphic  transformation  of  coordinates  7/,  :  — < •  5?"  with  ( Th)(x )  i— ►  (£,z),  with 
the  state  partition  in  the  new  coordinates  £  €  3?n_r,  z  £  3T  and  inverse,  (Ti,)(x)  (£,  z),  so 

that  the  full  state  representation  in  the  new  coordinates  is 


t  =  F(t,z),  (2.20) 

1  [-4(£,z)  +  B(£,z)u],  (2.21) 

where 

-4(£,~)  =  a(x)\x=f[t:) 

B((,z)  =  p(x)  |B=f(tiS). 


0  Ir-i 
0  0 


Definition:  The  zero  dynamics  of  the  input-output  model  (2.3)-(2.4)  are  given  by  the 
autonomous  system, 

i  =  F(t,  0).  (2.22) 

And  the  system  is  locally  minimum  phase  if  the  the  diffeomorphic  transformation  7\  is 
defined  on  a  neighborhood  of  the  origin  and  (2.22)  is  asymptotically  stable  to  the  origin; 
(  =  0. 


2.3  PLF  Computations  for  Nonlinear  MIMO  Systems 

The  nonlinear  input/output  model  is 

x  =  f(x)  +  G(x)u  (2.23) 

y  =  h(x)  (2.24) 

where  x  6  5?n,  u,y  6  5im  and  /,<#  (resp.  y.)  for  i  —  1, . . .  ,m  are  smooth  vector  fields  defined 
on  (resp.  3?m).  For  notational  simplicity  we  write, 

G{x )  =  [yi(x),...,ym(x)]. 

In  a  process  similar  to  the  previous  section  PLF  is  determined  by  transformation  of  the 
system  model  (2.23)-(2.24)  to  a  normal  form  in  which  we  identify  a  certain  ( mxm )  decoupling 
matrix  which  is  locally  nonsingular. 

The  process  begins  with  the  computation  of  an  appropriate  generalization  of  the  MIMO 
system  relative  degrees.  Let 

r,  :=  min{fc  =  1,2, :  Lgj(Lkf~l(hi))  7^  0,  for  some  j  =  1, .  .  . ,  m},  (2.25) 

the  ith  characteristic  number  [Fre75] .  Each  r*  is  then  the  minimal  relative  degree  of  the  set 
of  m  individual  output  responses  y*  obtained  from  each  input  uj  for  j  =  1, . . .  ,m.  Let 

'  <*!(*)  ' 

.  am(z)  . 


,4(x)  = 
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Nonlinear  PTant  Model 


_  _  _  _  _  _  _  I 


Figure  2.2:  Partial  Linearizing  Feedback  via  Nonlinear  Inverse  Model 


Figure  2.3:  Partial  Feedback  Linearization  and  Zero  Dynamics 


TSI-89-12-12-WB 


12 


where  a^x)  =  Lrf(h{)  and 


/5h(.t)  •••  0m{x) 


B(x)  = 


L  Pml{x)  0mm(x) 

where  f3ij{x)  =  Lgi(Lrf~^  (hi)).  Then  the  desired  normal  form  coordinates  are  2  E  #P  where 
r  =  j  r;,  which  are  obtained  as 

/  *1  \ 


z  = 


~2 


(2.26) 


V  / 

with  each  z;  G  9?r'  for  i  =  1, . . . ,  m  in  the  form, 


/  ^  \ 

Lf(h-i) 


V  Lr}-\h,)  ) 


(2.27) 


Proposition:  Given  the  system  (2.23)-(2.24),  there  exists  a  diffeonvorphic  transformation 
( T)(x )  =  (z,0  to  normal  form  coordinates, 


i  =  Az  +  E{A(z)  +  B(z)u} 
*  =  F{sti) 


(2.28) 

(2.29) 


where  A  =  diag{.di, . . . ,  i4m}  with 


A,  = 


0  lrx- 1 

0  0 


of  dimension  r,  x  r;  for  each  i  =  1, . . . ,  m  and  E  is  an  r  x  m  matrix  with  elements  given  as, 


(l,  if  i  =  r,-  and  j  =  i 

^ij  ~  \  0,  otherwise 

Then  the  PLF  control , 

u  —  -B~l(x){A(x)  -  i'}, 

(2.30) 

renders  the  v  t— ►  y  input/output  response  in  linear  form, 

z  =  Az  +  Ev, 

(2.31) 

II 

0 

u 

(2.32) 

Definition:  The  system  (2.23)— (2.24)  (output  constrained)  zero  dynamics  are  given  by 

Z  =  F(  0,0-  (2-33) 
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Definition:  We  say  the  system  is  locally  minimum  phase  if  the  zero  dynamics  are  asymp¬ 
totically  stable  to  the  origin  £  =  0. 

Remark:  In  general  the  computation  of  the  normal  form  with  (2.29)  independent  of  u  is 
difficult  and  not  required  to  establish  the  minimum  phase  property.  Instead,  note  that  the 
system  zero  dynamics  are  just  the  dynamics  of  (2.23)— (2.24 )  constrained  to  the  manifold 
Mh  C  5?"  of  dimension  n  —  m  given  by, 

Mh  =  {.T  G  R"  :  h{x)  =  0}. 

Proposition:  The  zero  dynamics  are  asymptotically  stable  if  and  only  if  the  system, 

x  -  f(x)  -  G(x)B~1(x)A(x),  z(0)  G  Mh 

is  asymptotically  stable  to  origin.  We  note  that  Mh  is  an  integral  manifold  {or  (2.23)-(2.24). 

2.4  Partial  Linearization  and  Variable  Structure  Control  Systems 

The  theory  of  Variable  Structure  (VS)  systems  addresses  the  design  of  control  laws  which  are 
discontinuous  functions  of  the  system  state.  VS  control  offers  practical  solutions  for  systems 
employing  actuators  which  can  be  efficiently  operated  in  bang-bang  and  other  discontinuous 
modes.  Our  interest  in  VS  control  for  rapid  slewing  of  multibody  systems  arises  from  the 
following  observations: 

1.  Design  methods  for  VS  control  for  output,  regulation  have  been  shown  to  effect  an 
implicit  partial  feedback  linearization.  The  implicit  partial  linearization  is  achieved 
through  the  use  of  VS  control  requires  only  output  feedback. 

2.  Direct  implementations  of  VS  control  laws  can  attain  a  level  of  robustness  to  plant 
model  assumptions  which  is  difficult  to  achieve  with  smooth  control.  Moreover,  ro¬ 
bustness  is  achieved  without  overly  conservative  restrictions  on  performance.  The 
limitations  of  robustness  for  VS  designs  are  clearly  related  to  the  minimum  phase 
conditions  for  PLF  design. 

3.  The  use  of  discontinuous  control  and  the  connections  between  VS  and  PLF  designs 
suggest  several  alternatives  for  integration  of  various  types  of  actuators,  including  both 
switched  and  continuous  modes  of  operation.  Integration  of  several  control  actuation 
systems  will  be  required  for  implementation  of  rapid  slewing  requirements  for  several 
candidate  large  space  platforms. 

The  general  theory  of  VS  control  design  is  well  known  and  we  will  not  attempt  to  present 
a  complete  description  of  its  scope.  For  details  see  the  survey  [DZM88].  However,  to  focus 
attention  on  the  concepts  we  seek  to  exploit  we  start  with  a  brief  description  of  the  basic 
ideas. 
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VS  control  systems  utilize  high  speed  switching  control  to  drive  the  system  trajectories 
toward  a  specified  manifold  called  the  switching  surface.  Given  the  nonlinear  system  (2.23), 
the  VS  control  laws  are  of  discontinuous  type; 


u{  = 


( 


for  s,(a;)  >  0 
for  s;(x)  <  0 


(2.34) 


with  s,(x)  —  0  smooth  switching  surfaces  chosen  in  the  state  space  for  each  i  =  1, . . .  ,m 
The  design  approach  which  is  preferred  is  based  on  the  introduction  of  sliding  modes. 


Definition  (sliding  modes):  A  manifold,  A4,,  consisting  of  the  intersection  of  p  <  m 
switching  surfaces,  s,(x)  =  0,  with  the  property  that  s,i,  <  0  for  each  i  —  l,...,p  in 
the  neighborhood  of  almost  every  point  in  M,  is  called  a  sliding  manifold.  Under  these 
conditions  any  trajectory  of  the  system  which  enters  M.,  remains  confined  to  the  manifold 
for  a  finite  length  of  time.  We  call  the  motion  on  M.,  a  sliding  mode. 

VS  design  methods  involve  a  two  step  process:  1)  design  the  switching  surface  so  that  once 
sliding  is  achieved  the  natural  sliding  mode  achieves  design  objectives  such  as  regulation, 
stabilization,  etc.,  and  2)  design  of  discontinuous  control  laws  which  achieve  sliding  on  desired 
regions  of  the  switching  surfaces.  The  method  of  equivalent  control  is  a  popular  approach 
for  designing  the  switching  surface  to  achieve  desired  sliding  mode  dynamics. 

Given  the  system  (2.23)  and  a  manifold,  M,  —  {x  €  :  s(x)  =  0},  with  s  :  3?"  — »  3Jm 

then  sliding  is  characterized  by  satisfaction  of  the  constraint  equations, 

s(x)  =  0,  i(x)  =  0  (2.35) 

over  the  finite  time  interval,  t\  >  t  >  t2  where  s(x(G))  =  0.  Note  that  a  sliding  mode  is  an 
instance  of  an  integral  manifold  for  the  closed  loop  system.  The  equivalent  control,  ueq,  is 

the  control  required  to  maintain  the  system  trajectory  within  the  manifold  M,  and  is  given 

by  the  condition, 

i  =  Vs(x)  x  =  Vs(x){/(x)  +  G(x)ue,}  =  0  (2.36) 

where  Vs(x)  =  ds/dx.  Under  the  assumption  that  det{Vs(x)G'(i)}  ^  0  for  x  G  Ais  we 
have, 

«e,  =  -[Vs(x)G(x)}-'Vs(x)f(x),  (2.3T) 

and  the  motion  in  sliding  is  given  by, 

x  =  {/  -  G(x)[Vs(x)G(x)]_1  Vs(x)}f(x),  s(x(G))  =  0.  (2.38) 

Connections  between  the  design  of  VS  control  and  feedback  linearizing  control  have  re¬ 
ceived  considerable  attention  [FH87].  The  principal  focus  has  been  on  the  problem  of  synthe¬ 
sis  of  VS  designs  for  nonlinear  systems  of  the  form  (2.23)  using  (exact)  feedback  linearization. 
In  the  sequel  we  direct  attention  to  the  problem  of  output  regulation.  The  connection  we 
establish  with  PLF  design  also  illuminates  several  questions  relative  to  robustness  properties 
of  VS  designs  with  sliding  modes. 
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Design  objective — Output  Regulation.  Given  the  system  (2.23)-(2.24)  where  y  indi¬ 
cates  a  set  of  regulated  outputs,  the  control  problem  is  to  drive  the  outputs  asymptotically 
to  zero. 

Since  output  regulation  problem  seeks  to  enforce  the  set  of  constraints 

h{(x)  =  0,  i  =  1, . . . ,  m, 

asymptotically  it  seems  reasonable  that  VS  design  could  be  employed.  However,  the  naive 
choice  3i(x)  =  hi(x)  leads  to  the  complication  that  in  general — in  fact,  most  often— [V/j(.t)G(x)] 
is  singular  for  almost  all  x  €  The  approach  suggested  in  [KK89]  is  to  design  sliding  mode 
via  the  choice  of  switching  surfaces  relative  to  the  normal  form  coordinates  for  ( 2.23 )—( 2.24 ) 
as  given  by  (2.28 )— ( 2.29). 

Proposition:  Given  the  system  ( 2.23 )— ( 2.24 )  obtain  the  diffeomorphic  transformation 
given  by  (2.26),  (2.27)  to  the  form  (2.28)-(2.29).  The  selection  of  switching  surface, 

s{z)  =  Kz  (2.39) 

with  K  an  m  x  777  constant  matrix,  solves  the  output  regulation  problem  if  sliding  can  be 
achieved.  In  sliding  the  equivalent  control  is, 

utq  =  -B{x)~1KAz  -  B~l(z)A(z)  (2.40) 

and  the  sliding  dynamics  are  given  by  the  r  linear  equations, 

:  =  [Ir-  EK]Az ,  A'r(0)  =  0  (2.41) 

Proof:  The  proof  is  given  in  detail  in  [KK89]  together  with  a  method  for  stabilization. 

By  defintion,  the  transformation  (T)(x)  *-+  (z,£)  is  invertible;  ( T)(z,( )  >->■  x,  and  the 
switching  surface  s(z)  =  0  can  be  reflected  to  the  original  coordinates  s(x)  =  0.  The 
dynamics  in  the  x  coordinates  are  best  understood  in  terms  of  the  geometry.  Define  the 
m-dimensional  manifolds  M,  =  (x  €  Sfn  :  $(x)  =  0}  and  Mh  =  {x  G  3?"  :  h(x)  =  0}.  The 
n  -  r  dimensional  manifold  Mz  =  {x  €  3?"  :  x  =  !T(0,£)}  is  contained  in  both  M,  and  Mh- 

Assume  that  in  some  neighborhood  V  €  3?n  a  sliding  mode  exists  on  V,  =  Pfl  M,  which 
is  assumed  nonempty.  Suppose  that  Vz  =  V,  ft  Mz  is  nonempty.  Let  a  denote  a  bounded, 
stable  attractor  of  the  zero  dynamics  contained  in  Vz.  Assume  that  all  trajectories  in  V. 
converge  to  a.  Then  if  the  initial  state  is  sufficiently  close  to  V,  the  trajectory  will  eventually 
reach  V ,  and  sliding  will  occur.  Clearly  the  stability  of  the  attractor  is  critical  to  the  stability 
of  the  overall  design.  In  the  sequel,  we  establish  conditions  for  output  regulation  of  multibody 
systems  which  guarantee  that  the  zero  dynamics  have  well  defined  local  equilibria  so  that 
linear  stability  analysis  of  the  zero  dynamics  is  appropriate.  We  remark  that  the  problem  of 
establishing  estimates  for  the  domain  of  attraction  in  the  zero  dynamics  is  an  open  question. 
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2.5  Robust  Stabilization  of  Nonlinear  Systems 

For  practical  implementation  of  PLF  compensation  various  researchers  have  focused  atten¬ 
tion  on  conditions  which  guarantee  robust  stabilization  of  the  nonlinear  system  (2.3)-(2.4) 
with  feedback  transformation  of  the  form  (2.11)  by  introduction  of  linear  feedback  v  =  Kz. 
A  brief  survey  of  the  wide  range  of  methods  which  have  been  proposed  is  given  in  [BBKA88]. 
In  this  section  we  focus  attention  on  a  ubiquitous  assumption  in  most  of  the  work  on  robust 
stabilization  of  nonlinear  systems  based  on  PLF  compensation. 

Consider  the  usual  case  for  engineering  design  where  the  open  loop  system  dynamics  for 
( 2.23 )-( 2.24)  is  given  by  a  nominal  model  of  the  form 

x  =  f°{x)  +  G°{x)u  (2.42) 

y  =  h°(x),  (2.43) 

where  f° , g°  are  C°°  vector  fields  for  i  =  l,...,m,  defined  on  a  manifold  At  G  9i",  with 
f°( 0)  =  0.  We  assume  the  (true)  system  response  can  be  modeled  via  a  perturbation  of  the 
vector  fields; 

f  =  f°  +  A/,  G  =  G°  +  AG 

with  A/,  A <7,  each  C 00  defined  on  M  and  A/(0)  =  0.  In  [AB88]  a  detailed  analysis  is  given 
leading  to  sufficient  conditions  on  A/,  A G  which — together  with  the  assumption  that.  (2.42) 
is  feedback  linearizable — guarantees  that  (2.23)  is  also  linearizable.  The  conditions  given  are 
less  restrictive  then  the  usual  structure  matching  conditions  [AB87],  The  structure  matching 
conditions  also  play  a  role  in  establishing  conditions  for  robust  stabilization  and  we  repeat 
them  for  convenience. 

The  structure  matching  conditions.  Under  the  assumption  that  the  perturbation  vec¬ 
tor  fields  satisfy, 

A/,  Agi  e  A  (2.44) 

where  A  =  5p{<7i, .  .  .  ,<7m}  then  any  such  model  (2.23)  is  exactly  linearizable — in  fact,  by  the 
same  diffeomorphtc  transformations.  The  above  conditions  are  equivalent  to  the  statement 
that,  the  perturbations  to  the  vector  fields  can  be  factored  as: 

A/(x)  =  G°(x)df(x),  (2.45) 

A  G(x)  =  G°{x)Dg(x).  (2.46) 

The  importance  of  the  structure  matching  conditions  in  establishing  robust  stability  is 
that  under  these  conditions  the  model  uncertainty — after  application  of  PLF  compensation — 
can  be  equivalently  represented  by  a  perturbation  (or  disturbance)  at  the  compensated  sys¬ 
tem  inputs.  This  facilitates  the  design  of  disturbance  rejection  techniques  using  either  ex¬ 
plicitly  nonlinear  control  designs  such  as  Gutman  [GL76]  or  linear  control  design  such  as  in 
[Kra87j.  To  see  this  we  summarize  the  construction  under  the  assumption  that  h(x)  =  h°(x). 
Substitute  the  PLF  compensation  (2.30)  obtained  for  the  nominal  model  (2.42)-(2.43)  into 
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the  model  for  the  true  system  (2.23)  to  obtain, 

x  =  f°  +  Af  +  {G°  +  AG)B~l[v  -  A] 

=  f°  +  G°Adf  +  G°{I  +  Dg)B~l[v  -  A] 

=  f0  +  G°B-l[v-A}  +  Go[df  +  DgB(v  -A)}.  (2.47) 

The  model  after  nominal  PLF  compensation  can  (by  the  above  assumptions)  be  transformed 
to  r-coordinates  defined  in  (2.26)-(2.27)  to  obtain, 


z  =  Az  +  £{v  -f  7;}  (2.48) 

i  =  (2-49) 


where 


[<!/(*)  +  D,(i)B(»r,{t,  -  ^(x)}]|i=r.,(!i 

[d,(x)  -  D,(*)B(*)-1>1(i)]|„t.,(s)  +  [ca(x)B(x)-'»] 


(2.50) 

(2.51) 


Thus  it  is  clear  that  robust  stabilization  must  address  the  disturbance  rejection  of  the  class 
of  input  disturbances  dv(t)  which  bound  the  model  error;  ||t/||  <  \\dv\\.  An  effective  design 
approach,  in  the  case  when  the  nominal  model  is  exactly  feedback  linearizable,  is  given  by 
Spong  and  Vidyasagar  [SV87].  There  approach  utilizes  an  L stabilization  criterion  and 
obtains  a  linear,  time-invariant  feedback  control  for  the  v-input. 


Remarks:  The  design  methods  for  robust  stabilization  of  nonlinear  systems  available  in  the 
literature  are  almost  exclusively  based  on  the  structure  matching  conditions.  By  reflecting 
the  model  uncertainty  to  the  system  inputs  (after  nonlinear  feedback  compensation)  they 
can  employ  either: 

1.  linear  compensator  design  which  seeks  to  reduce  loop  gains  consistent  with  bounded 
model  uncertainty,  or 

2.  nonlinear  switching  mode  compensator  design  which  seeks  to  over-bound  input  distur¬ 
bances  by  high  gain  implementations  using  fast  switching  control. 

Both  methods  result  in  essentially  conservative  designs  since  the  worst  case  bounds  on  the 
input  disturbances  must  be  assumed. 

In  [Akh89]  the  basis  for  robust  stabilization  of  nonlinear  systems  is  considered  further 
and  new  results  are  obtained  for  the  case  of  parametric  model  uncertainty.  The  new  con¬ 
trol  laws  obtained  in  [Akh89]  employ  basic  constructions  of  adaptive  control  in  the  context 
of  feedback  linearization.  The  results  show  that  robust  stabilization  can  be  obtained  under 
much  less  restrictive  conditions  than  the  structure  matching  conditions.  Significantly,  feed¬ 
back  linearization  plays  the  role  of  enforcing  linearity  for  subsequent  control  loop  designs. 
To  the  extent  that  this  can  be  achieved  in  practical  applications  it  can  enhance  reliability 
and  repeatability  thus  achieving  improved  performance  prediction — an  important  feature 
for  space-based  systems.  Standard  constructs  in  adaptive  control  can  then  be  applied  to 
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enhance  the  robustness  of  feedback  linearization  with  model  uncertainty.  We  emphasize 
that  the  constructions  described  below  are  new  and  offer  stability  results  for  the  nonlinear 
system  under  very  general  assumptions.  In  the  next  few  paragraphs  we  briefly  review  some 
significant  aspects  of  these  results. 

Robust  Stabilization  of  Nonlinear  Systems  by  Adaptive  Methods.  Again,  starting 
with  the  system  model  in  the  form  (2.23)  we  assume  the  model  uncertainty  can  be  represented 
by  parametric  dependence  of  the  vector  fields  so  that  the  model  has  the  form, 

i  =  f(x,0)  +  G(x,6)u,  (2.52) 

with  9  a  p-vector  of  unknown  parameters.  We  assume  that  for  every  9  £  Be,  a  closed, 
compact  neighborhood  of  the  nominal  parameter  0o ,  /  and  gt,  for  i  =  1  ,...,m,  are  G°° 
vector  fields  and  f(0,9)  =  0  for  9  £  Be.  The  nominal  design  model  is  characterized  by  the 
set  of  nominal  parameters  and  we  take  /(.r)  =  f(x,90),  G(x)  =  G(x,90). 

The  following  assumptions  are  used  in  [Akh89]  to  establish  a  robust  stabilizing  controller 
for  the  nonlinear  system. 

Assumption  1:  The  nominal  system  is  exactly  feedback  linearizable.1 

It  is  important  to  note  that  assumption  1  is  applied  only  at  the  nominal  plant  model 
with  fixed  and  known  parameters  90.  This  is  in  contrast  to  the  structure  matching  conditions 
which  are  almost  universally  assumed  in  the  current  literature. 

Assumption  2:  For  any  x  £  U  C  and  9  £  Be, 

A gl(x,9)  £  Sp  {gi{x,90), . .  .,gm(x,90)} 

for  i  =  1, . . . ,  m.  This  assumption  implies  that  there  exists  anmxm  matrix  valued  function, 
D(x,6),  with  smooth  elements  such  that  AG  =  G°D. 

Assumption  3:  Either  there  exists  an  m  x  m,  strictly  positive  definite  matrix,  K  such 
that, 

Vx  £  U,  VO  £  Be,  0  <  D(x,0)  <  K , 
or  there  exist  a  K ,  negative  definite,  such  that, 

Vx  £  [/,  V0  £  Be,  K  <  D{x,0)  <  0. 

This  assumption  (in  various  forms)  is  typical  in  adaptive  control  stability  analysis  and 
design.  It  says  that  the  “sign”  of  this  term  must  be  definite  and  known  a  priori. 

The  design  approach  is  natural  and  begins  with  the  transformation 

-  =  T(x,e0), 


‘The  extension  of  the  results  described  below  to  the  case  of  stabilization  by  PLF  is  in  progress. 
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to  normal  form  and  choice  of  feedback  linearizing  compensation  (2.30).  The  v  control  is 
chosen  in  two  parts.  First,  for  the  nominal  design  and  performance  objectives  we  find 
v  =  Fz  where  A  +  EF  in  (2.31)  is  a  stable  matrix.  In  this  case,  there  exists  a  unique, 
positive  definite,  symmetric  solution  to  the  Lyapunov  equation, 

(.4  +  EF)rP  +  P(A  +  EF)  =  -I. 

The  design  for  the  nominal  model  is  now  modified  by  the  introduction  of  adaptation. 
The  control  law  obtained  is  described  by  the  following  equations; 


u  = 

-B  ^z)^®)  -  v}, 

(2.53) 

z  = 

T(x,0o ), 

(2.54) 

V  = 

Fz  +  C(zJ)6, 

(2.55) 

9  = 

—  T-1Ct(z,6)EtPz , 

(2.56) 

II 

-pB(z)GT(z,9o)Pz0T. 

(2.57) 

In  these  control  laws  the  p  x  p  matrix  T  >  0  is  the  “adaptation  gain”  which  is  chosen  so  that 
the  Lyapunov  function, 

v  =  zTPz  +  eTre 

is  positive  definite  for  all  x  €  U  and  p  is  a  positive  scalar.  The  result  established  in  [Akh89] 
is  the  following. 

Theorem:  The  adaptive  control  with  feedback  linearizing  transformation  is  asymptotically 
stable  to  the  origin  x  =  0  and  9  — >  9,  asymptotically,  if  for  x  £  U  and  9  £  Be  there  exists 
ci  >  1  and 

2||  TT(x)P<t(xJ)  -  wTDy  ||f>||!  ||  <  c2||T(i)||!, 

where 

<P(z,0)  =  [&f{z,0)  +  &G{z,6){A(z)  +  B-'Fz}} 

and 

y  =  GT(z,  90)PT z. 


Proof:  [Akh89]. 

2.6  Partial  Linearization  for  Lagrangian  Systems 

Despite  the  apparent  simplicity  of  determining  the  zero  dynamics  from  the  normal  form  as 
above  it  is,  however,  quite  complex  to  compute  the  complete  transformation  leading  to  the 
full  state  normal  form  as  given  above.  One  approach  (if  possible)  is  to  obtain  the  full  state 
exact  linearizing  transformation  via  the  procedure  given  by  Hunt,  Su,  and  Meyer  [HSM83] 
which  requires  the  solution  of  a  set  of  simultaneous  partial  differential  equations.  However, 
in  many  special  cases  the  zero  dynamics  as  well  as  the  required  transformations  for  partial 
linearizing  control  can  be  obtained  more  directly.  In  the  sequel  we  discuss  the  required 
constructions  for  Lagrangian  systems. 
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Consider  the  case  of  a  square  Lagrangian  system  with  inputs  r  £  and  outputs  y  £  5Rm. 
Suppose  that  the  n  generalized  coordinates  can  be  partitioned  into  components  qx  £  !Km  and 
92  £  so  that  the  equations  of  motion  take  the  form 


d 

dL 

dL 

dt 

dqi 

dqx  ~  T' 

(2.58) 

d 

dL 

dL 

dt 

dq2 

-  «-  =  o, 

dq2 

(2.59) 

y  -- 

=  h{qltq2). 

(2.60) 

Assume  that  the  origin  is  an  equilibrium  point  with  r  =  0,  h( 0,  0)  =  0,  and  that  the  Jacobian 
dh/dqi  is  nonsingular  on  some  neighborhood  of  the  origin.  Furthermore,  we  assume  that 
the  Lagrangian  is  a  positive  definite  quadratic  form  in  the  generalized  velocities.  Then  the 
input-output  map  (r  — *  y )  has  relative  degree  2  (locally),  a  PLF  control  exists  and  the 
zero  dynamics  may  be  computed  by  a  relatively  simple  coordinate  transformation  applied  to 
(2.58)— (2.60). 

In  order  to  demonstrate  these  properties  we  introduce  a  change  of  coordinates  (91,72)  1 ■ 

(y,u)  via  the  relations 

y  =  M  91,92),  u  =  92  (2.61) 

Note  that  the  assumption  detdh/dqx  ^  0  at  the  origin  assures  that  this  is  a  valid  local 
coordinates  transformation  and  the  inverse  relations  can  be  given  as 


q\  =  g{y,u),  q2  =  u. 


(2.62) 


Since  any  “point"  transformation  preserves  the  Lagrangian  structure  of  the  equations,  in  the 
new  coordinates  we  can  write  the  variational  problem  in  the  form 


d  dL  dL  - 

dt  dy  dy  yT' 

d  dL  dL  - 

dt  dii  du  u 


where 


and 


L{y,y,  u,u)  =  £(91,91,92-92)! 


qi=g(y,u) 


Ty  ~  [dy 

Equations  (2.63)-(2.64)  reduce  to  the  form 


dg 


i  t 


,  ru  = 


dg 

du 


jyy  +  Nii  +  I\y(y,y,u,u)  =  tyr , 
NTy  +  Juu  +  A'u(y,  y,  u,u)  =  fur. 


(2.63) 

(2.64) 

(2.65) 


(2.66) 


(2.67) 

(2.68) 
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Let  us  define  the  partitioned  matrices  Sy,  Su,  4?,,,  via  the  relations 

(*v\  =  \  Jy  *  V1  (  A'v  \ 

Vs-/  NT  UJ 

and 

)  _  jy  N  r1/  fy  \ 

)  nt  Ju  {  fu  )  ■ 

Note  the  choice  of  control 

reduces  (2.67)— (2.68)  to 

y  =  v  (2.72) 

NT(y.  u)y  +  Ju(y,u)u  4-  I<(y,y,u,  it)  =  rti(j/,u)$"1(i/.u){Sy(j/,u)  +  r}  (2.73) 

where  we  have  explicitly  displayed  the  dependence  of  the  model  parameters  on  the  generalized 
coordinates.  Equation  (2.72)  provides  the  linearized  input-output  dynamics  and  the  zero 
dynamics  are  obtained  from  (2.73)  upon  setting  y(t)  =  0,  which  implies  y  =  0,  y  =  0,  and 
v  —  0.  Thus,  we  obtain  the  zero  dynamics  in  the  form 

ju(0,  u)u  +  A'(0,  0,  u,  u)  -  fu(0,  u)$y  *(0,  u)Sy(0,  ti)  =  0  (2-74) 

which  represents  an  autonomous  nonlinear  dynamical  system  in  the  state  coordinates  i t,u. 
We  say  the  system  is  locally  minimum  phase  if  the  origin  in  the  (u,  »)  coordinates  is  a 
stable  equilibrium  for  (2.74).  If  the  system  is  minimum  phase  then  selecting  the  control  v 
to  stabilize  the  origin  of  (2.72)  guarantees  stability  of  the  origin  of  the  dynamical  system 
(2.58)— ( 2.60).  Thus  the  computational  complexity  of  obtaining  the  zero  dynamics  depends 
on  the  complexity  of  the  required  inverse  relation  g  in  (2.62). 


21 

(2.69) 

(2.70) 

(2.71) 
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3  Reduced  Order  Modeling  of  Multibody  Systems  and  Nonlin¬ 
ear  Control  Design 

The  approach  described  in  [BBKA88]  to  obtain  finite  dimensional  models  for  multibody  sys¬ 
tems  with  elastic  interactions  utilizes  Finite  Element  Methods  (FEM)  based  on  the  method 
of  collocation  by  splines.  In  this  section  we  focus  attention  on  a  system  model  for  attitude  re¬ 
orientation  including  a  primary  set  of  rigid  bodies  whose  attitude  relative  to  an  inertial  frame 
is  given  by  a  set  of  attitude  parameters  xT  =  The  system  model  includes  spatially 

distributed  elastic  interactions  with  FEM  deformation  coordinates,  x,  =  {rjT )T  £  9?2n'. 

After  some  manipulation,  the  general  class  of  multibody  systems  we  will  consider  can  be 
written  in  the  generic  form, 

Mr(xr,X,)xr  +  Nx,  +  Kr(zr,  Xr,Xs,Xs)  -  Tb  (3.1) 

Ma{xr,x,)x,  +  NTxr  +  A\(£r,xr,x„x,)  =  Gf,,  (3.2) 

where  rb  £  9?"'’  are  independent  torques  applied  to  rigid  bodies  and  f,  £  9?2fc  are  generalized 
forces  acting  on  the  flexible  structure.  By  virtue  of  the  collocation  FEM  model  coordinates 
we  note  that  the  influence  of  structural  control  enters  through  the  2 n,  x  2k  matrix, 


G  = 


Gf  0 
0  Gm 


where  G/,G'm  are  n,  x  k  with  elements  given  by, 

[G/hi  = 

v  ~ 

z i  control  moment  is  applied, 


1  if  at  z  —  Zi  control  force  fj  is  applied, 
0  else 


if  at 
0  else 


[C-lii 

As  a  result  there  exists  a  n,  x  n,  nonsingular  permutation  Va  such  that 

which  defines  a  change  of  basis  for  the  structural  deformation  coordinates, 


(3-3) 

(3.4) 


hk 

0 


with  x,ti  £  $?2fc,  xs<2  €  9?2(n,_fc).  The  model  can  then  be  decomposed  as, 

Mr{xr,x,)xr  +  NiX$ti  +  iV2*,,2  +  Kr(Xr,Xr,X,,Xs)  =  Tb, 
Mf'HXfi  -(-  Ms\2Xs2  +  Aj  .Tr  -b  A  ji  (  XT ,  XT ,  X  j ,  X  s  )  —  fs , 
+  Mt,22Xt,2  +  Nj  Xr  +  K.,2(xr,X  r,Xs,X,)  =  0. 

Let  m.  =  nb  +  2 k  and  introduce  the  definitions, 

/ 

Tb 


(3.5) 

(3.6) 

(3.7) 


(;:)■ 


6  :=  Xf.2 
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Then  we  can  write  the  model  (3.5)— (3.7)  in  the  simplified  form, 

+  A/126  +  —  /.  (3.8) 

hhiii  +  M 22^2  +  A2(£,£)  =  0.  (3.9) 

To  illustrate  the  practical  simplicity  of  decoupling  control  computations  and  the  fun¬ 
damental  design  issues  we  focus  attention  on  a  ?n-vector  of  independent  system  outputs 
consisting  of  linear  combination  of  deformation  coordinates, 

V  =  [CuC2]  V  (3.10) 

where  C  —  [C'i,  C'2]  is  constant,  m  x  n  matrix  with  n  =  nb  -f  2na.  Without  loss  of  generality 
we  assume  the  m  x  m  matrix  C\  is  nonsingular.  To  identify  the  specific  form  of  decoupling 
control  for  this  special  case  we  rewrite  (3.8)— (3.9)  in  ^-coordinates  (3.10); 

xy  +  (M12  —  il/uCj  lC2)i2  +  hy{y, y> £21^2)  =  />  (3.11) 

MnCx  ly  +  (A/22  —  M2\CX  1C2)^2  +  A2(t/, y, £2* £2)  =  0.  (3.12) 

It  is  then  straightforward  to  identify  the  required  decoupling  control; 

/  —  +  (T/i2  —  M 1 1 C\  lC2)^2  +  M 1 1  C\  1  u .  (3.13) 

Applying  the  decoupling  control  (3.13)  obtains  the  first  m  degrees  of  freedom  (3.11)  in 
decoupled,  linearized  coordinates  with  synthetic  control  in  ‘acceleration  coordinates’; 

V  =  u.  (3.14) 

The  remaining  n  —  m  degrees  of  freedom  (3.12)  comprise  the  dynamics  which  are  then 

decoupled  from  the  output.  We  note  that  the  £2  coordinates  will  be  driven  by  the  synthetic 

control  u.  However,  for  most  practical  designs  the  principal  consideration  is  for  the  stability 
of  the  zero-output,  constrained  dynamics  in  £2;  i.e. ,  the  system  zero  dynamics.  Here  the  zero 
dynamics  are  readily  identified  from  (3.12)  with  y  =  y  =  y  =  0; 

(M22  —  1C'2)^2  +  A2(0, 0,  Z21Z2)  =  0.  (3.15) 

The  form  of  (3.15)  is  significant  since  PFL/decoupling  is  feasible  only  if  the  zero  dynamics 
are  stable  in  an  appropriate  sense.  In  this  study  we  take  rapid  slewing  to  be  a  ‘rest-to-rest’ 
maneuver.  At  the  end  of  the  maneuver  we  will  require  precision  pointing;  i.e.,  structural 
alignment.  Thus  we  focus  attention  on  asymptotic  stability  of  the  structural  alignment 
state;  i.e.,  £2  =  0,£2  =  0  in  (3.15).  As  seen  from  (3.15),  a  general  choice  of  outputs  as 
linear  combination  of  the  deformation  coordinates  together  with  the  primary  body  attitude 
coordinates  may  lead  to  nonlinear,  possibly  unstable  zero  dynamics. 

For  various  practical  considerations  we  focus  attention  on  implementation  of  decoupling 
control  for  the  special  choice  of  outputs  collocated  with  the  control  forces ;  i.e., 

y  =  {inh,GT}  (  )  =  I7-0!  (  f2  )  • 


(3.16) 
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This  choice  of  outputs  for  decoupling  will  obtain  the  zero  dynamics  in  the  form, 

M22Z2  +  -^2(0, 0,  £2,  £2)  =  0.  (3-17) 

More  importantly,  our  assumptions  on  elastic  potential  energy  for  the  structural  deformation 
[BBKA88]  focuses  attention  on  small  amplitude  deformation  dynamics  and  we  have, 

Ks(xr,xr,x„  x,)  =  Bcx,  +  Kcx,  +  D(xr,xr,xs,  xs)  (3.18) 

with  D  :  5fn  — ►  5UV'  with  D  =  0  for  xT  =  0  (i.e.,  when  the  rigid  body  is  at  rest),  and 
Bc ,  Kc  constant,  positive  semidefinite  matrices  (see  Section  5.1).  Under  these  assumptions  it 
is  easy  to  show  that  (3.17)  is  a  set  of  n  —  m  linear,  second  order  equations  for  the  structural 
deformation  dynamics  constrained  at  the  physical  output  locations  given  by  (3.10);  i.e., 
the  structural  dynamics  with  certain  combinations  of  localized  pinned  and/or  cantilevered 
supports.  Thus  stability  properties  of  the  zero  dynamics  for  the  case  of  outputs  collocated 
with  the  control  forces  follow  from  the  natural  structural  properties  of  stiffness  and  damping. 

3.1  PFL  with  Reduced  Order  Models 

In  this  section  we  will  outline  general  considerations  for  implementation  of  PFL  via  nonlinear 
partial  state  feedback.  Considerations  for  rapid  slewing  of  flexible  space  structures  suggest, 
that  the  relative  stiffness  of  the  structure  by  comparison  with  the  severity  of  the  maneuver 
dynamics  will  dominate  the  complexity  of  the  process  model  for  control  design.  In  addition  to 
nonlinear  couplings  the  distributed  parameter  nature  of  the  structural  deformation  suggests 
that  time  scaling  via  singular  perturbation  methods  be  applied.  In  the  previous  section  we 
described  nonlinear  controller  realization  which  achieves  exact  decoupling  using  full  state 
feedback.  We  remark  that  considerations  for  good  mechanical  design  suggest  that,  space 
systems  designed  for  rapid  slewing  maneuvers  will  be  structurally  stiff  [Le87].  Considerable 
attention  has  been  given  to  the  use  of  composite  material  for  enhanced  damping  and  stiffness 
properties  of  such  structures  [RFL87] .  However,  a  primary  system  cost  for  any  space  system  is 
mass  and  the  resulting  constraints  on  system  mass  will  ultimately  limit  structural  stiffness. 
Our  models  have  been  configured  to  represent  generic  considerations  for  rapid  slewing  of 
space  structures  but  have  been  scaled  to  focus  dynamic  effects  expected  to  dominate  behavior 
of  optical  systems  such  as  the  SBL. 

In  previous  studies  where  the  singular  perturbation  approach  was  used  [Dwy88,  KO88], 
a  scaling  parameter  e  was  introduced  as  a  measure  of  the  system  compliance;  for  example, 
Dwyer  [Dwy88j  considers  e  =  ^  while  dimensionless  scaling  of  the  form  e  =  jjjj-  [KO88],  is 
more  often  chosen.  In  both  cases,  scaling  is  chosen  relative  to  uq,  the  lowest  modal  frequency 
due  to  the  structural  flexure,  and  u>0,  the  natural  frequency  of  the  rigid  system.  For  relatively 
stiff  systems  the  scaling  leads  to  e  <  1.  Using  standard  singular  perturbation  approach,  the 
dynamics  of  the  system  can  then  be  approximated  by  the  dynamics  of  the  reduced  slow 
subsystem  (i.e.,  rigid  body  dynamics). 

We  have  considered  a  slightly  different,  time  scale  separation  which  identifies  a  reduced 
order  model;  viz.,  a  quasi-rigid  model.  The  dynamics  of  the  quasi-rigid  model  do  not  coincide 
necessarily  (as  in  previous  time  scale  separations)  with  the  rigid  body  dynamics.  The  slow 
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modes  will  consist  of  the  rigid  body  modes  together  with  an  arbitrary  number,  p,  of  low 
frequency  flexure  modes.  The  fast  modes  will  consist  of  the  remaining  (higher  frequency) 
flexure  modes.  Thus  we  can  study  the  tradeoff  of  model  order  reduction  vs.  decoupling 
performance  directly. 

We  assume  that,  any  structure  designed  for  repeated  rapid  slewing  will  be  relatively  stiff 
and  we  investigate  the  obvious  time  scale  decomposition  of  (3.1 )— ( 3.2)  as  follows.  We  assume 
that  Ms  is  positive  definite  and  Kc ,  Bc  in  (3.18)  are  each  positive  semidefinite.  To  obtain 
appropriate  time  scaling  we  transform  (3.1 ) — (3.2)  to  modal  coordinates  as  follows.  Under 
the  above  assumptions  there  exists  a  nonsingular  matrix  P  such  that: 

PtAI,P  =  In,  and  PTI<CP  =  K 

where  K  =  diag(a;i, . . .  )  with  u>i  <  u>2  <  . . .  <  uin,-  We  also  assume  P  is  such  that 

PT  BCP  =  B 

where  B  is  a  diagonal  matrix. 

Note  that  xj  —  PTxa  transforms  (31),  (3.2)  from  displacement  coordinates  x,  into  modal 
coordinates  X/; 

M(xr,Xf)xr  +  NT'if  +  A'r(xr,  xT,Xf,if)  =  rb  (3.19) 

'if  +  Nxr  +  Bci}  +  Kcxf  +  D{xT,xr,Xf,Xf)  =  Gf,  (3.20) 

where 

D  =  PTD(xr,ir,PTx„PTi,)P, 

N  =  NP,  G  =  PtG,  M{xr,Xf)  =  Mr(xr,  PTxs) 

and 

7?r(xr,xr,x/,x/)  =  A'r(xr,xr,  PTx„  PTxs).  (3.21) 

In  the  sequel  the  overbars  will  be  omitted  and  we  assume  the  model  is  given  in  modal 
coordinates,  (3. 19)— (3.20). 

Time  scale  decomposition  is  obtained  by  identifying  a  separation  between  the  first  p  low 
order  modes,  u>J,...,u>p  from  the  remaining  higher  modal  frequencies  u>p+1 , . . . ,  .  The 

separation  above  induces  a  decomposition  of  the  state  Xf  as: 

x,=  X\  ,  x)  G  Rp ,  x)  €  Rn-~p, 

.  xt 

with  the  corresponding  decomposition  for  the  matrices  N,  G,  Kc,  Bc  and  D  being 


II 

r,  with 

No  €  Rp, 

Ni  G  Rn’~p 

(3.22) 

II 

c5 

C) 

-■s 

r,  with 

Go  e  Rpx2k , 

G\  G  A(n'“p)xm 

(3.23) 

*■  -  ( *.■  i ) 

|  ,  with 

Ko  €  Rpxp, 

A’i  G  A^n'-p^x^n’_p^ 

(3-24) 

|  ,  with 

Bo  G  Rpxp , 

B\  G  /?(n--p>x(n’_pl 

(3.25) 
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Xr  ,  Xr ,  Xy,  Xy  )  ^  ^  (  Xr  ,  Xr ,  Xy ,  Xy  ) , 

with  Dx  :  5?n  -  3?p,  £>2  :  5T  - 


(3.26) 


With  this  decomposition,  (3.20)  becomes: 

Xy  4-  N0Xr  4-  Boi°f  +  Kox°f  4-  Di(.T°,.Ty)  =  Gif,  (3.27) 

x/1  4-  N\XT  +  B\Xfl  4-  K\X^  +  Z?2 ( ^ / >  x})  =  Gif,.  (3.28) 

By  assumption  the  modal  frequencies  u>*+1, . . . ,  u/^  are  of  the  same  (large)  order  of  mag¬ 
nitude  and  we  can  express  them  as  multiples  of  1/e  with  e  =  Thus  the  stiffness, 

A’i  =  diag(ujp+1 ,  o,’p+2, .  . .),  can  be  scaled  as: 

A'i  =  ^A'io,  (3.29) 

and  the  “fast  time”  state,  x y,  will  satisfy, 

Xy  =  ez  where  z  ~  0(1).  (3.30) 

Substituting  (3.27),  (3.28),  (3.29)  and  (3.30)  into  (3.19),  (3.20)  yields  the  standard  form  for 
singular  perturbation  analysis; 

Mr(xr,  x°,  ez)xr  4-  N^x°f  +  eNjz  +  A'r(xr,  xr,  x°,  x°,  ez,  ez)  =  rh  (3.31) 

x°f  4-  N0xr  +  B0x°f  +  A'0Xy  4  Dx{xr,  x°,  ez)  =  Gif,  (3.32) 

e[z  +  B\z]  +  Nixr  +  Di(xr,  x°j,  ez)  +  Iiioz  =  Gif,.  (3.33) 

The  reduced  state  consists  of  the  rigid  state,  xr,  and  a  p- dimensional  part  of  the  flexible 
state  for  the  elastic  structural  deformation. 


3.1.1  Quasi-Rigid  Model 

The  singular  perturbation  approach  models  the  time  responses  as  decomposed  into  slow  time 
or  quasi-steady  state  and  boundary  layer  terms.  Neglecting  the  boundary  layer  amounts  to 
letting  e  — *  0  and  the  resulting  slow  subsystem  is: 


M0(xr,x°f)xr  +  Nq  x°f  +  A'r°(xr,xr,x°,x°) 

=  Xj,, 

(3.34) 

x°  4-  N0xr  +  B0x°f  +  A'0Xy  4-  D°(xr,  xr,  x°,  x°) 

II 

Q 

> 

(3.35) 

NiXr  +  A'io~  +  £>°(xr,Xr,Xy,Xy) 

=  Gif,, 

(3.36) 

where 

M0(xr,x°f,ez)  =  M(xr,x°f,0) 

and 

A’°(xr,  xr,  x”,  Xy,  ez,  ez)  =  A'r(xr,  xr,  x°f,  x" 

,0,0). 
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The  dimension  of  the  state  space  is  thus  reduced  from  2 (nr  -f  na)  to  2(nr  +p)  and  we  identify 
the  quasi-rigid  model  in  the  form 


where 


Mo{xr,x0s ) 


T° 

Xf 


+  C0(Xr,  Xr,  X°fl  X°)  =  E 


M  0(Xr,X0f) 
C0(xr,  Xr,  X  j,  X  j  ) 

E 


Mo{xr,x°f)  jV0r 

No  Ip 

K?(xr,xr,x°/}x°f) 

B0x°f  4-  K0x°f  4-  D°(xr,x7.,x°,x°) 

In  0  \ 

o  6'i  / 


(3.37) 


(3.38) 

(3.39) 

(3.40) 


By  assumption  A'10  is  nonsingular  and  the  n,  -  p  algebraic  equations  (3.36)  yields  the 
following  expression  for  z  in  terms  of  xr,xr,  x°,x°  and  rj,: 


-K^lDliXr.Xr^X^i))  +  NXXr  ~  G2fa], 

-A'fo1  [D°2(xr,xr,x0f.i°f)  +  Ni(M0  -  NoNlr'in  -  A'rn(Tr, ir, i°) 
+Nq  (Boi°f  +  K0x°f  +  D°(xr,xr,  x°,  x°f)  -  GJa )}  -  G2fa ] . 


(3.41) 


Equation  (3.41)  defines  a  2(nr  +  p)-dimensional  manifold  M0  in  the  2(nr  +  na)  dimensional 
state  space  called  the  slow  manifold. 

The  quasi-rigid  model  (3.37)  approximates  the  slow  response  as  a  quasi-steady-state  for 
the  full  model  (3.31 )— (3.33).  The  difference  between  the  response  of  the  quasi-rigid  model 
and  that  of  the  full  model  (3.31 )— ( 3 . 33 )  is  given  by  a  the  boundary-layer  system  which  is 
obtained  as  follows.  Let  zj  be  the  fast  time  scale  part,  of  z  (i.e.,  in  the  stretched  time  scale 
r  =  ~7j).  Then  the  fast  time  scale  system  is; 


dr2 


+  B\o~~  +  hio~f  =  G2f! 


(3.42) 


where  the  Bl0  is  obtained  via  the  appropriate  scaling  of  the  damping;  B i  =  yfeB i0. 

The  trajectories  of  the  full  system  (3.31 )— (3.33)  can  be  approximated  by  examining  the 
solutions  of  the  quasi-rigid  model  (3.37).  The  approximation  is  0(e)  under  the  assumption 
that  the  fast  subsystem  (3.42)  is  asymptotically  stable.  For  e  small  (i.e.,  sufficient,  separation 
of  adjacent  modal  frequencies)  stability  of  the  boundary  layer  must  be  obtained  by  inherent 
damping  or  by  the  introduction  of  fast  time  scale  (wide  bandwidth)  control  //.  This  anal¬ 
ysis  suggests  the  alternatives  for  material  damping  vs.  active  control  will  be  significant  in 
achieving  robust  high  performance  nonlinear  decoupling  control.  The  role  of  stabilization 
of  the  boundary  layer  is  in  improving  the  approximation  offered  by  model  order  reduction. 
The  forces  required  are  relatively  weak  but  of  wide  bandwidth.  They  can,  in  principle,  be 
generated  by  internal  material  properties  deliberately  introduced  through  the  use  of  passive 
damping  or  by  active  control  using  imbedded  actuators. 
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3.1.2  The  reduced  flexible  model 


In  analyzing  the  qualitative  robustness  of  nonlinear  decoupling  and  PFL  we  may  require 
refinements  of  the  quasi-steady  state  analysis  described  above.  In  this  section,  we  describe 
alternatives  for  refined  reduced  order  modeling  by  application  of  the  method  of  integral  man¬ 
ifolds.  Following  [KK086],  we  note  that  the  solutions  x(t,  e),  z(t,  e)  of  ( 3 . 30 )— ( 3 . 32 )  consist 
of  a  fast  boundary  layer  and  a  slow  quasi-steacly-state.  The  boundary  layer  is  significant 
only  in  z(t,e)  since  x(t,e)  is  predominantly  slow  and  its  boundary  layer  is  no  larger  than 
0(e).  As  we  have  seen  in  the  previous  section,  (3.41),  obtained  for  e  =  0,  defines  a  2(n,  -f-p)- 
dimensional  manifold,  A/o  C  9?n.  For  nonzero  €,  we  define  a  2(n,  +  p)-dimensional  manifold, 
M(  C  depending  on  the  scalar  parameter  t  by 


M. 


6  5?"  : 


h(x,x,rb,  fa,e),  z  =  h{x,x,Tb,ft 


<e)J  - 


(3.43) 


where  it  is  assumed  that,  h  is  continuously  differentiable  sufficiently  many  times  in  all  of  its 
arguments.  The  manifold  Mt  is  said  to  be  an  integral  manifold  for  the  system  (3.44)-(3.30) 
if.  for  given  initial  conditions  r(fl),  ,r(fl),  r(0),z(0)  in  A/e,  then  the  trajectories  x(<),  z(t.)  are 
in  A/e  for  all  t  >  0.  Following  [KK086],  it  follows  that,  the  existence  of  a  stable  equilibrium 
manifold  A/0  of  ( 3.34 )—( 3.36)  for  e  =  0  (since  A'10  is  nonsingular),  implies  the  existence  of 
an  integral  manifold  A/t  of  (3.31)-  (3.33).  When  the  fast,  dynamics  (3.42)  are  asymptotically 
stable,  then  if  z(0),  i(0)  are  not  initially  in  A/«,  z  will  converge  to  Mt  after  the  decay  of  a  fast 
transient;  Zf.  Thus  the  response  x(t)  of  the  full  system  (3.31 )— (3.33)  will  rapidly  approach 
M(  and  then  flow  along  A/e.  As  e  — >■  0,  the  manifold  A/e  converges  to  M0. 

By  definition,  the  function  h  defining  the  integral  manifold  in  (3.43),  satisfies  (3.33) — the 
manifold  condition; 

e  [h  +  Bih ]  4-  Nxxr  +  h\0h  +  D2(xT,xr,x0f,x°f.eh,eh)  =  G^f,-  (3.44) 


Solving  for  h  in  (3.44),  then  by  replacing  z  by  h  and  z  by  h  in  (3.31)— (3.32),  we  obtain  the 
reduced  flexible  model  [SKK87]: 

MT{xr,eh)xr  +  NqX0,  -(-  eNfh  +  Kr{xr,xr,eh,eh)  =  rb  (3.45) 

ic°f  +  N0xr  +  B0x°f  +  K0x°f  +  Di(x°f,x°f,eh,eh)  =  G\f,-  (3.46) 


We  remark  that  although  this  model  has  the  same  dimension  as  the  quasi-rigid  model  (3.37) 
it  is  not  an  approximation  to  the  singularly  perturbed  model  ( 3.31 )— (3.33 );  instead,  it.  repre¬ 
sents  the  exact  system  (3.1 )— ( 3.2)  restricted  to  the  manifold  A/e.  The  reduced  flexible  model 
(3.45 )— (3.46)  therefore  models  the  flexible  system  more  accurately  than  does  the  quasi-rigid 
model  (3.37).  Unfortunately,  to  construct  the  reduced  flexible  model  ( 3.45)— (3.46),  one  needs 
(in  general)  to  solve  a  partial  differential  equation  for  h.  The  approach  presented  in  [SKK87], 
is  to  approximate  the  manifold  Mt  and  the  reduced  flexible  model  up  to  any  order  in  e.  The 
first  step  consists  of  expanding  h(x,x,e)  in  a  power  series; 


h(x,  x,  e)  =  h 0(r,  ,t)  +  ehflx,  x)  +  .  .  . 


(3.47) 
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and  substituting  into  the  manifold  condition  (3.33)  to  obtain  the  coefficients.  Then  the 
control  is  synthesized  via  an  expansion; 

r5(x,i,e)  =  r°(x,x)  +  er^(x,i)  +  . . .  (3.48) 

f,{x.x,e)  =  fS0(x,i)  +  ef3l{x,x)  +  ...  (3.49) 

which  highlights  the  importance  of  the  relatively  weak  forces  in  the  refined  analysis.  To 

compute  the  hi  s,  we  equate  terms  in  like  powers  of  e  in  the  manifold  condition  (3.33)  from 
which  we  obtain 

h0  -  -Kw[D2(xr,Xr,x0f,i°f)  +  Ni£r  -  G2f,„]  (3.50) 

h i  =  -A io*  [TMxr,  xr,  ho,  ho)  +  B\ho  4-  ho  —  G^/,,  ]•  (3.51) 


As  expected,  h0  has  the  same  expression  as  (3.41).  We  can  also  write: 

M(x,  eh )  =  A/0(x)  +  cA/i(x,  h0)  +  . . .  (3.52) 

A'r(x,  x,  eh,  eh)  =  A‘r°(x, x)  +  eA',!(x,  x,  h0,  h0)  +  ■  ■  ■  (3.53) 

Substituting  (3.47)-(3.49),  (3.52),  (3.53)  in  (3.45)— ( 3.46)  obtains  a  reduced  flexible  model 
up  to  any  desired  degree  of  accuracy  in  e.  Let  the  m-dimensional  control  vector  /  =  (rr,  fj  )T 
and  recall  that  x  =  [x^,x”T]T.  The  reduced  flexible  model  (3.45)— (3.46)  can  then  be  written 

[M0(xr)xr  +  Co(xriXr)  ~  E f0] 

+  e  i ( xr ,  h-Q )xr  +  C i(xr,  x r i  h0,  h0 ,  h0)  -  Efi ]  +  0(e2)  =  0  (3.54) 

where  M0(xr),  C0(xr,xr)  and  E  are  as  defined  in  (3.38)~(3.40).  The  order  1  correction 
terms  Ai\  and  C\  are 


Mx{xr,h0)  = 
C\(xr,  xT,  ho,  ho,  ho) 


MX(Xr,ho)  0 

0  0 


(3.55) 

(3.56) 


Kl(xr,  xr,h o,  h0)  +  Nf  h0 

D^ho) 

Thus  neglecting  terms  of  order  e2  and  higher  in  (3.54),  we  obtain  the  First-Order  Corrected 
Model  (FOCM)-, 

[Mo  +  eM\){xr,  h0)xr  +  [C0  +  eCi)(xr,  xr,  h0,  h0,  h0)  =  E(f0  +  efx),  (3.57) 
an  0(e)  improved  approximation  to  the  reduced  flexible  model. 


3.1.3  Feedback  linearization  and  decoupling  of  the  quasi-rigid  model 

Recall  the  quasi-rigid  model  (3.37): 

Mo{x)x  +  C0(x,  x)  =  Efo  (3.58) 

where  x  =  [xJ,x°fT]T  and  /0  =  [r°  T ,  ffn}T ■  Depending  on  the  number,  2k,  of  structural 
controllers  available  and  the  number  p  of  flexible  states  available  for  feedback,  the  controls, 
and  f3,  can  be  synthesized  using  feedback  linearization  to  completely,  or  partially,  decouple 
the  quasi-rigid  model  (3.58)  from  the  deformation  z.  We  consider  three  cases  as  follows. 
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Case  1  (m  <  p):  An  important  question  which  we  first  consider  concerning  the  ap¬ 
proximating  model  (3.58)  is  whether  it  is  exactly  feedback  linearizable?  That  is,  can  the 
linearization  be  performed  by  identifying  a  critical  set  of  outputs  for  which  the  input/output 
relation  is  invertible?  If  such  a  set  of  outputs  can  be  found  then  no  further  consideration  for 
decoupling  and/or  stability  of  system  zero  dynamics  is  required.  To  answer  this  question, 
let  us  write  the  model  in  state  space  form  with  Xi  =  x  =  [xf,x°T]T  £  Rn’+P, 

(  ^  )  =  d(xux2)  +  G(xi,x2)f  (3.59) 


where 


d(x  i,jc2)  = 


z2 

-Mo\xx)C0{xi,x2) 


)' 


GM  =  Mo:(ta)E. 


(3.60) 


Conditions  for  exact  linearization  are  well  known  but  can  be  tedious  to  check.  The 
required  computations  are  straightforward  and  can  be  performed  using  a  symbolic  algebra 
system.  We  utilized  CONDENS,  a  symbolic  manipulation  package  for  analysis  and  design 
of  nonlinear  control  systems  using  geometric  methods,  to  check  the  feedback  linearizability 
conditions  [?].  Checking  the  conditions  for  each  p  =  0, 1,  2, . . .,  we  find  that  the  first,  condition 
(i.e.,  controlability)  is  generically  satisfied;  i.e., 

dim{S,  [/,  G\,ad2fG, . . . ,  ad)p+lg}  =  2p  +  2. 


However,  as  we  suspected,  if  m  <  p  the  second  condition  of  linearizability — the  involutivity 
condition  [Isi85j — is  not  satisfied.  This  means  that  the  distribution 

Sp  {g,[f,G)...,ad?g} 

is  not  involutive  and  hence,  that  the  system  (3.58)  is  not  exactly  feedback  linearizable  in 
the  case  m  <  p.  This  is  the  generic  case  for  control  of  distributed  parameter  systems  with  a 
finite  number  of  localized  controls. 

To  consider  PFL  we  identify  a  set  of  m  system  outputs  for  the  inverse  dynamics  com¬ 
putations.  A  natural  choice  which  guarantees  the  minimum  phase  property  is  to  identify 
collocated  outputs  as: 

y  =  Etx. 

In  this  case  the  decoupled,  quasi-steady  state  zero  dynamics  include  the  dynamics  of  the 
last  p  -  m  modes  of  the  structure  with  additional  boundary  conditions  arising  from  y  =  0. 
After  the  decay  of  the  fast  transients,  the  uncontrolled  deformation  z  approaches  the  slow 
manifold  M0  defined  by  (3.41)  with  the  controls  and  f,  replaced  by  the  resulting  linearizing 
input-output  controls. 

Case  2  (m  =  p):  In  this  case,  the  reduced  order  model  (3.58)  is  exactly  feedback  lineariz- 
able  with  the  linearizing  control 


/o  =  E  1  (C0(,r,j)  +  Afo(a-)v)  • 


(3.61) 
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where  v  =  (v‘[,vJ)T  E  Rm .  The  effective  quasi-steady  state  response  from  the  synthetic 
control  v  to  output  y  is  given  by: 

x  =  v.  (3.62) 

As  in  the  previous  case,  the  deformation  2  approaches  the  slow  manifold  Mo  defined  by  (3.41) 
with  r°  and  /°,  the  components  of  /0  are  as  in  (3.61). 


Case  3  (m  >  p):  If  the  number  of  dominant  flexible  modes  p  we  choose  to  incorporate 
in  the  reduced-order  model  is  smaller  than  the  number  of  available  structural  controls,  then 
p  of  these  controllers  can  be  used  for  decoupling  and  linearization  of  the  reduced  order 
model  (3.58),  while  the  remaining  m  —  p  controllers  can  be  used  to  shape  the  slow  manifold, 
Mo,  as  well  as  to  compensate  for  the  structural  elastic  response,  2,  within  reasonable  limits. 
The  ideas  is  due  to  Dwyer  [Dwy88]  who  studied  the  method  of  deformation  shaping  in  the 
particular  case  p  =  0;  i.e.,  when  the  reduced  order  model  (quasi-rigid  model)  coincides  with 
the  rigid  body  dynamics  and  the  “fast”  subsystem  coincides  with  the  dynamics  of  the  whole 
flexible  appendage. 

In  the  case  m  >  p,  the  structural  control  forces,  /,,  can  be  decomposed  as 


/.= 


fs,\ 

/j,  2 


f 3,1  e  Rp ,  f 3,2  e  R2k~p. 


The  p  x  m  full  rank  matrix  Gi  in  the  slow  subsystem  (3.34)— (3.36)  can  be  put  in  the 
form  [Gu  0]  using  elementary  permutations,  where  Gn  is  a  p  x  p  and  nonsingular  matrix. 
Similarly,  the  (2 k  —  p)  x  m  matrix  G 2  takes  the  form  [0  ,  G22]  with  the  (2k  —  p)  x  (m.  —  p) 
matrix  G22  having  full  rank  (m  —  p).  The  reduced  order  model  (3.58)  can  then  be  written  as 


Mo(x)x  +  C0(x,x)  =  EifZ 
where  /»  =  (r6°  r,/5°/)r  E  Rn‘+P  and 


Ei  = 


//  0 
V  0  Gn 


(3.63) 


(3.64) 


The  model  (3.63)  is  exactly  feedback  linearizable  with  the  linearizing  feedback  as  in  (3.61) 
but  with  E 1  replacing  E.  The  algebraic  constraint  (3.41)  defining  the  slow  manifold  M0 
becomes  in  this  case 


2  =  -K1-01(D2(xr,ir,x°f,i°f)  -f  +  G22f,„_r],  (3.65) 

where  the  acceleration  xr  can  be  replaced  by  the  synthetic  control  tq.  Contrary  to  the 
previous  cases,  the  deformation  2  can  be  controlled  or  compensated  (completely  in  the  case 
m  =  q)  by  the  structural  controllers  /,m  which  will  require  high  bandwidth,  but  “low- 
authority”  (i.e.  relatively  weak)  actuators  for  implementation. 
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3.1.4  Partial  feedback  linearization  of  the  Reduced  Flexible  Model 

In  some  applications,  and  depending  upon  the  choice  of  the  scaling,  e  may  not  be  “small” 
enough  relative  to  the  decoupling  requirement.  In  this  case,  it  is  necessary  to  replace  the 
quasi-rigid  model  (3.58)  by  the  more  accurate  FOCM  (3.57)  for  control  design.  The  appli¬ 
cation  of  feedback  linearization  to  the  FOCM  introduces  a  correction  term  in  the  controls 
which  improves  the  decoupling  (and  thus  linearization)  of  the  system  by  an  order  of  e.  The 
corrected  controls  are  obtained  via  an  expansion; 

f  =  f°  +  cf\  (3.66) 

where  f°  is  the  linearizing  torque  for  the  quasi-rigid  (zero-order  approximation)  model. 

Design  of  the  corrective  control  fl :  The  corrective  control  fl  is  designed  to  improve 
the  decoupling  by  annihilating  the  0(e)  terms  in  (3.57),  i.e., 

[A4i(.t,  h0)x  +  Ci(x,  x,  h0,  h0,  h0)  -  F/1]  =0.  (3.67) 

The  construction  will  depend,  as  in  the  previous  section,  on  the  values  of  p  and  m.  In  the 
case  m  <  p  the  corrective  control,  /*,  can  be  constructed  to  annihilate  part  of  (3.67).  In  the 
case  m  =  p,  (3.67)  can  be  exactly  annihilated  by: 

fv  =  E~l  ( Ci(x ,  x,  h0,  ho,  ho)  +  Mi(x,  h0)vj  .  (3.68) 

If  we  apply  the  control,  /  =  f°  +  e/1,  to  the  reduced  flexible  model  (3.54),  where  r”,  rj ,  /° 
and  fl  are  as  designed  above,  we  obtain: 

i  =  v  +  0(e2)  (3.69) 

which  linearizes  the  reduced  flexible  system  up  to  order  e2. 

3.2  Conclusions  and  Directions 

The  results  described  in  this  section  represent  a  preliminary  analysis  of  the  problem  of  PLF 
compensation  for  multibody  system  with  flexible  interactions.  (A  preliminary  version  of  the 
analysis  was  presented  at  the  1989  American  Control  Conference  in  Pittsburgh,  PA  [ABB89]. 
Our  interest  in  this  phase  of  the  research  was  primarily  motivated  by  concerns  for  the  com¬ 
plexity  of  implementation  of  PLF  compensation  required  for  high  order  FEM  models  of  the 
elastic  response  of  a  multibody  system.  The  primary  objective  was  to  demonstrate  rapid 
slewing  and  precision  pointing  of  the  system  attitude  as  defined  relative  to  a  body-fixed  frame 
attached  to  the  primary  system  body.  The  analytical  approach  based  on  integral  manifold 
methods  offers  a  refinement  of  the  standard  singular  perturbation  approach  which  suggests 
that  higher  order  corrections  applied  to  the  rigid  body  torques  can  improve  performance  of 
PLF  compensation  for  relatively  low  order  models  of  structural  flexure.  Simulation  studies 
performed  during  this  study  did  not  reveal  substantial  improvement  from  such  corrections. 

In  response  to  the  observed  responses  of  simulation  studies  we  can  draw  several  conclu¬ 
sions.  The  SBL  system  model  developed  in  this  study  consists  of  a  multibody  system  with 
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primary  system  outputs  to  be  regulated  consisting  of  rigid  body  positions  and  parasitics 
arising  from  low  mass  structural  interactions.  Primary  system  controls  are  collocated  with 
the  position  coordinates  for  regulation.  Singular  perturbation  analysis  reveals  that  the  sys¬ 
tem  flexure  dynamics  in  such  cases  are  only  weakly  controllable  in  the  sense  that  control  of 
the  fast  time  scale  enters  only  through  the  slow  time  system.  As  is  well  known,  such  system 
will  require  additional  control  authority  to  affect  control  of  the  fast  time  scale  response. 
For  the  problem  of  rapid  slewing  and  precision  pointing  this  means  that  active  structural 
control  which  interacts  with  the  low  mass  density,  structural  components  may  be  important 
in  achieving  enhanced  robustness  of  the  system  slewing  and  alignment  responses.  This  ob¬ 
servation  is  potentially  important  given  the  increasing  reliance  on  passive  damping  methods 
for  space  structure  applications  [RR87],  Unfortunately  a  comprehensive  study  will  require 
some  analysis  of  emerging  technology  for  imbedded  actuation  and  sensing  of  flexible  struc¬ 
ture  response.  This  was  viewed  outside  the  scope  of  the  second  year  efFort  on  the  subject 
contract. 
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4  Some  Design  Approaches  for  Combined  Rapid  Slewing  and 
Precision  Pointing 

Consider  the  general  form  of  the  equations  of  motion  for  a  multibody  system  undergoing 
attitude  reorientation.  The  model  includes  n  degrees  of  freedom  (including  elastic  degrees 
of  freedom  retained  from  a  FEM  expansion)  with  position  coordinates  q  £  5?n  and  velocities 
p  €  :Rn.  The  model  has  the  form, 

q  =  f( q )p  (4.1) 

M(q)p  +  B{q,p)p  +  K(q,p)q  =  Gf  (4.2) 


where  M,  B ,  A"  are  n  x  n  symmetric,  matrix- valued  functions  and  G  is  an  n  x  m  matrix. 
The  m- vector  /  includes  generalized  forces  (and  torques)  acting  on  the  multibody  system. 
For  simplicity  we  assume  the  position  coordinates  are  ordered  as:  q  =  ,  u7']7'  where 

7  G  9i3  parametrizes  the  rotation  of  the  body  fixed  frame  of  the  system  principal  body. 
Also,  p  —  [u>r,  iiT]T ,  where  u>  is  the  body  angular  rates.  The  general  model  includes  nonlinear 
kinematics  (4.1)  which  is  assumed  to  have  the  form, 

r(7)  o 

0  In  — 3  . 

where  the  3  x  3  matrix  is  defined  previously  depending  on  the  parametrization  of  50(3). 

The  general  problem  of  control  of  system  pointing  and  multibody  alignment  is  addressed 
by  establishing  a  set.  of  m  principal  system  outputs  taken  as  an  7n-vector  of  position  coordi¬ 
nates. 

V  =  Cq.  (4.3) 


f(?) 


Claim  1:  The  decoupling/PLF  control 


/  =  (CTM~lG)~1v  -  (CT M~lG)~1CT D  -  B)p  -  Kq } 


(4.4) 


with 


D(  7)  = 


1TvI3  03tn_3 
— 3,3  On— 3,n  — 3 

renders  the  input/output  response  from  v  — *  y  as  a  linear,  2m-dimensional  system  of  the 
form, 

y  =  v. 


Proof:  To  obtain  (4.4)  we  proceed  by  direct  computation  of  the  normal  form  coordinates  for 
the  general  model  (4.1)— (4.2)  with  outputs  (4.3).  First,  identify  a  change  of  state  coordinates, 

m(<?, ?)-(--, a 
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=  Om,n  C 
.  Z3 

with  Z\,Z->  chosen  so  that  T  has  rank  2??.  Thus  a  direct  computation  yields, 


Cq 

(4.6) 

C'q  =  Ct(q)p 

(4.7) 

Z\q  +  Z\V(~i)p 

(4.8) 

an  invertible  transformation. 

Before  we  obtain  the  system  model  in  the  (r,£)  coordinates  we  express  the  accelerations 
in  a  convenient  form.  From  (4.1 )— (4.2)  we  have, 

d  d 

i  =  jti  =  +  r(<7)p.  (4.9) 

Note  that  the  n-vector  valued  function, 


'  J{rhW  0 
0  0 


r(7)  0 

0  In  — 3 


(4.10) 

(4.11) 


A  principal  advantage  of  the  use  of  the  Gibbs  vector  parametrization  for  attitude  position 
is  the  simple  algebraic  form  of  the  kinematic  relation.  Appendix  A  establishes  the  following 
based  on  the  Gibbs  vector  parametrization, 


d_ 

d~t 


{r(7Wr(7)w  =  7ru>r(7)u;. 


Using  this  relation  we  obtain, 


<7 


7rwr(7)o; 

'  n7)  0 

0 

0  In  — 3  . 

M(q)-'{Bp+I<q-Gf} 


t(y)M(q)-'{\HHq)D(q)- B(q,p)]p- K(q,p)q  +  C.f} 


where, 

iTvi3  0 

0  0 


The  model  can  be  expressed  in  normal  form  coordinates  as, 

“1  =  “2 1 


22  =  + 

i  = 


(4.12) 

(4.13) 


(4.14) 

(4.15) 

(4.16) 
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From  (4.5),  (4.13)  we  have, 

A<1  >P)  =  f(7)M(^)_1{(il/(^)Z)((7)  -  B{q,p))p- K(q,p)q}  (4.17) 

B(q,p)  =  T(i)M(q)~1G.  (4.18) 

The  decoupling  control  is  then, 

/  =  ^(P>9)_1{i'  ~  -4(/>,?)}. 


and  (4.4)  follows. 


□ 


Claim  2:  Given  the  model  (4.1 )— (4.3)  we  can  find  ( n  —  m )  x  n  valued  functions  Z\(q),  Z2(q) 
such  that 


and 


is  nonsingular. 


Zlt(q)M(q)-lCT  =  0, 
Z2t(q)M(q)-lCT  =  0, 


C 


T{q)  ~ 


0 

.  Zi(q) 


o  ' 
c 

Zi{q)  . 


(4.19) 

(4.20) 


Claim  3:  In  the  case  that  outputs  for  regulation,  linearization,  and  decoupling  defined  by 
(4.3)  are  collocated  with  the  generalized  control  forces;  i.e. , 

span{G}  =  span{C,r} 


then  the  (output-constrained)  zero  dynamics  are  given  as, 


£=  F(*.*)la=0 


where  F(z,£)  =  $(p,q)  with  (p,q)  such  that 


and 


nmq)  (;)«(;) 


${p,q)  =  Zx{q)q  +  Zl(q)t(q)p  +  Z2{q)t{-f)p 

+  Z2(q)r(q)M-1{[M(q)D(q)-  B(q,p))p+  K(p,q)q}. 


(4.21) 


(4.22) 


(4.23) 
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Proof:  From  (4.13)  and  (4.5)  we  obtain  £  =  Zxq+  Z2q.  Direct  evaluation  of  (4.16)  obtains, 

£  =  Zl(q)q  +  Zt(q)q  +  Z2(q)t('y)p  (4.24) 

+  Z2(q)f(‘f)M(q)-1{[M{q)D{q)  -  B{q,p))p  -  I<{q,p)q  +  Gf}.  (4.25) 

Then  by  Claim  2  and  the  assumption  we  have, 

Z2(q)t(1)M(q)-lG  =  0. 

a 


4.1  Practical  Implementation  of  Time-optimal  Slewing  for  an  Inertial  Load 

The  advantage  of  the  PLF  control  which  implements  the  system  inverse  is  the  slewing  control 
can  now  be  implemented  in  the  new  coordinates.  A  switching  control  law  is  designed  for  the 
independent  decoupled  axis  motions  in  acceleration  coordinates.  The  decoupled  single  axis 
response  obtained  after  PLF  compensation  is: 

0{t)  =  a{t)  (4.26) 

where  ft  is  the  attitude  parameter  for  ith  axis  (subscripts  suppressed)  and  a  (the  acceler¬ 
ation)  is  the  control.  Requirements  for  rapid  slewing  suggest  a  minimum  time  criterion  is 
appropriate  for  the  design  of  control  laws. 

Necessary  conditions  for  a  control  law,  optimal  with  respect,  to  a  minimum  time  objective, 
can  be  obtained  directly  from  the  maximum  principle.  Under  magnitude  constraint  on  the 
available  acceleration, 

HOI  <  «n.M  (4.27) 

the  optimal  control  is  of  bang-bang  type  with  switching  strategy  described  by  a  switch¬ 

ing  surface  in  the  state  space  9?2  with  coordinates  x  =  [0,0]r.  For  the  double  integrator 
plant  the  optimal  control  involves  determining  a  single  switch  point.  Briefly,  the  optimal 
switching  surface  is  a  quadratic  curve  in  5?2  which  can  be  obtained  as  follows.  For  constant 
accelerations,  a(t)  =  amRX,  the  natural  trajectories  are, 

0(0  =  6+tw  t,  (4.28) 

ft(0  =  to  +  tit  +  (4.29) 

where  0(0)  =  £0  and  0(0)  =  0  are  initial  conditions.  The  optimal  trajectory  reaches  the  equi¬ 
librium  at  the  origin  in  minimum  time.  To  identify  the  required  switching  surface  choose  the 
trajectories  which  reach  the  origin  at  time  tf  by  solving  the  above  equations  simultaneously 
and  eliminating  dependence  on  tf,  to  obtain, 

tf  =  2<w&.  (4.30) 

It  is  easy  to  see  that  the  required  switching  surface  in  the  state  space  3?2  with  coordinates 
z  =  [0,  ft]T  is 


3(2)  =  f(z  2)  +  2amnxzi  =  0, 


(4.31) 
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with  /(~2)  =  ~2 1 ~2 1  and  the  time-optimal  control  is  bang-bang; 

a{t)  =  —  amnx  sgn{s(z)}. 


(4.32) 


However,  any  practical  implementation  of  the  time-optimal  control  law  will  involve  an  ap¬ 
proximation  of  the  ideal  discontinuous  actuation.  Most  realistic  actuators  have  slew  rate 
limits,  dead  zones,  delays,  or  other  dynamics  which  may  play  a  role  in  implementation  dy¬ 
namics.  A  simple  example  of  the  difficulty  of  implementing  time-optimal  control  can  be 
illustrated  if  we  replace  the  ideal  signum  function  with  a  saturating  element.  In  this  case 
it  can  be  shown  that  the  origin  is  not  an  asymptotically  stable  equilibrium.  To  see  this 
consider  a  Lyapunov  function, 


defined  on  5?2  with  coordinates  as  given  above.  For  the  practical  control  law, 


where 


q(0  = 

:  amaxsat{s(z)}, 

(4.33) 

, 

f 

for  t  >  l/tfi 

sat.(£)  :=  | 

9it, 

for  |£|  <  l/gx 

(4.34) 

{  -1, 

for  £  <  -l/gi 

where  gx  is  a  positive  constant  gain  representing  slew  rate  limiting  in  the  actuator. 

In  a  neighborhood  of  the  origin  such  that  jz2|  <  1  / gx  we  can  evaluate  the  time  derivative 
of  the  Lyapunov  function  along  trajectories  of  the  closed  loop  system; 


V'  =  zTb 

-  J 


O 

1  ■ 

'  o ' 

[0 

o 

i 

sat{2amBXr1  4-  z2\z2\} 


0  1 
-2amnx  |  ~2 1 


(4.35) 

(4.36) 


(4.37) 


The  V-function  is  positive  definite,  decreasent,  and  radially  unbounded.  In  the  neighborhood 
of  the  origin  we  have  V  <  0  for  \z2\  <  l/<7i-  Any  finite  precision  implementation  of  f{z2) 
will  ultimately  become  negligible  near  the  origin  and  the  closed  loop  system  response  will 
ultimately  oscillate  with  frequency  \/2a:mBX  and  amplitude  given  by  the  precision  of  /  near 
the  origin. 

The  above  limitations  are  simple  examples  of  the  general  tradeoffs  involved  in  design 
of  nonlinear  control  systems  where  dual  mode  operation  must  be  considered  to  resolve  the 
tradeoffs  between  requirements  for  static  accuracy  such  as  local  stability,  noise  rejection,  etc. 
and  dynamic  accuracy  such  as  speed  of  response.  These  issues  are  well  known  to  control 
engineers  but  do  not  receive  a  lot  of  attention  in  the  theoretical  literature.  Considerations 
for  dual  mode  operation  in  this  case  suggest  that  the  optimal  switching  function  be  replaced 
with  a  dual  mode  version  which  approximates  time  optimal  operation  for  large  signal  motion 
but  achieves  asymptotic  stability  of  the  origin  for  small  angle  displacements.  Additionally, 
considerations  for  small  signal  overshoot,  stability  margins,  etc.  can  be  included  without 
sacrifice  of  the  large  signal  response. 
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A  dual  mode  version  of  the  time-optimal  control  law  can  be  obtained  by  replacing  the 
nonlinear  function  /  with  a  peicewise  continuous  modification  with  finite  slope  near  z2  =  0. 
The  choice, 


f(=2)  ■= 


,2 

“2 

+  9A 
^  4  ’ 

04 

Al 

€M 

*2 

52^2, 

for  |z2|  <  g2/2 

,2 

“2 

£l 

4  ’ 

for  z2  <  -g2/ 2 

(4.38) 


obtains  asymptotic  stability  of  the  origin  for  any  g2  >  0  with  closed  loop  control  (4.33). 

The  resulting  practical  implementation  of  the  time-optimal  servo  control  law  is  shown 
in  Figure  4.1.  Here  0C  is  the  commanded  attitude  angle.  To  the  extent  that  practical 
limits  of  actuators  permit,  the  gains  g\  and  g2  should  be  chosen  to  satisfy  local  stability 
and  overshoot  as  well  as  large  angle  rapid  (time-optimal)  slewing.  These  tradeoffs  can  be 
reconciled  by  appealing  to  the  ideal  bang-bang  control  law  (4.32)  with  the  modified  /.  With 
ideal  switching  the  surface,  s(r)  =  0,  (4.31)  is  an  integral  manifold  for  the  closed  loop 
system.  From  the  theory  of  discontinuous  differential  equations  we  can  obtain  stability 
analysis  by  decomposition  of  the  closed  loop  system  dynamics  into  two  phases:  1)  reaching 
to  the  switching  surface  and  2)  ideal  sliding  on  the  switching  surface. 

Assume  the  closed  loop  system  achieves  ideal  sliding  condition  s(z)  =  0  for  t  >  t,.  Then 
the  sliding  dynamics  are  obtained  from  the  constraint  i(r)  =  0  for  t  >  ta; 


•s(“)  —  f'(:  z)z2  +  2amaxZi  —  0 


where  /'  =  df/dz2.  In  the  linear  region  of  the  peicewise  function  / — namely,  when  |z2|  < 
g2/ 2 — we  obtain, 

s(^)  =  92&  +  2amnxz2  =  0.  (4.39) 

The  equivalent  control  in  sliding  is 


2am(VX 

—  Z2. 

92 


(4.40) 


The  sliding  dynamics  on  the  manifold  M,  =  {z  £  :  s(z)  =  0}  are  then  obtained  by 

substitution  of  the  equivalent  control  in  the  closed  loop  equations.  In  the  linear  region  we 
obtain, 


2  ar 


•2  - 


92 


and  together  with  the  constraint  s(z)  =  0  the  ideal  sliding  dynamics  reduce  to 

2am.x 


9  2 


■~2 


(4.41) 

(4.42) 


(4.43) 


which  are  asymptotically  stable  to  z2  =  0  for  any  g2  >  0  The  local  time  constant  in  sliding 
is 


r  = 


92 


2«r 


(4.44) 
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To  choose  the  gain  g2  we  take  #x  as  the  desired  tolerance  for  large  angle  slewing;  i.e. ,  for 
any  initial  offsets  exceeding  0x  we  want  to  apply  large  angle,  time-optimal  slewing  control, 
but  for  any  offsets  smaller  than  0X  we  want  linear  mode  action.  Then  the  required  gain  is 
obtained  by  substituting  zi  =  9X  into  (4.31); 

92  =  (4.45) 

The  resulting  linear  mode  time  constant  is 

r  =  --  — ■ 

y^mnx 

Note  that  the  ideal  sliding  dynamics  retain  the  feature  of  no  overshoot  characteristic  of 
time-optimal  slewing  into  the  linear  mode  operation. 

To  implement  the  time-optimal  servo  we  need  to  understand  the  limits  of  saturation 
mode  operation  of  the  actuator  and  the  effect  of  available  torque.  From  the  theory  of  sliding 
mode  design  it.  remains  to  show  in  what  region  sliding  will  be  achieved.  A  straightforward 
application  of  the  reaching  analysis  of  [DZM88J  follows.  From  (4.31)  let  V(z)  =  |s2.  To 
achieve  sliding  we  need  to  show  that  ss  <  0.  Here 


i(z)  =  f'{z  2)z2  +  2amnxz1 


(4.46) 


where 


2z2,  for  z2  >52/2 
<  g2,  for  | z2 1  <  g2/2 

.  -2z2,  for  z2  <  -(72/2 


(4.47) 


A  necessary  and  sufficient  condition  for  sliding  to  be  achieved  globally  with  the  ideal  switch¬ 
ing  control  law 

<*(0  =  -amnxSgn{s(z)}, 


is 


^  ^niRx  • 


Within  the  linear  region  sliding  is  achieved  when 


max  • 


For  practical  implementation  of  time-optimal  slewing  the  choice  of  actuator  slew  rate  gx 
is  critical.  With  the  (practical)  closed  loop  control  (4.33)  the  dynamics  in  a  neighborhood 
of  the  origin  are  given  by 


0  1 
— 2am(>x<7i  ~otmaxgxg2 


(4.48) 


Sensitivity  of  the  stability  properties  can  be  determined  from  a  root,  locus  analysis  with 
respect  to  gx  (see  Figure  4.2).  The  characteristic  equation  of  the  closed  loop  system  is 


^  T  otm*xgxg2\  +  2anl(ix<7i  —  0 
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Figure  4.1:  Implemenation  of  Time-Optimal  Servo 


Figure  4.2:  Root  Locus  for  Linear  Mode  Sensitivity  of  Saturation  Mode  Slewing 
which  we  write  in  the  form 

A2  +  k(g2\  +  2am(VX)  =  0 

with  k  =  amKXgi .  To  approximate  ideal  sliding  in  the  linear  region  we  require 

8 

9\  >  ~ 

9  2 

which  guarantees  that  the  roots  are  real  and  stable.  As  gi  >0  is  increased  the  small 
signal  response  is  dominantly  a  single  time  constant  which  approaches  the  ideal  sliding 
mode  response  in  the  limit  as  gx  — *  oo  [SS83] . 

4.2  Implementation  of  Rapid  Slewing  with  Coordinated  Modes  of  Actuation 

In  this  section  we  describe  the  basis  for  design  of  coordinated  use  of  different,  modes  of  ac¬ 
tuation  for  rapid  slewing.  The  idea  is  to  utilize  the  smooth  (or  continuous)  mode  actuators 
(such  as  AMED’s)  to  implement  a  component  of  the  PLF  compensation.  This  implemen¬ 
tation  permits  a  simple  decomposition  of  the  slewing  control  into  individual  single  axis 
commands.  The  slewing  control  is  implemented  using  discrete  (discontinuous)  mode  actua¬ 
tion  (such  as  reaction  jets).  This  is  consistent  with  the  very  large  torque  requirements  for 
large  angle,  rapid  slewing  of  significant  system  inertias.  Bang-bang  mode  operation  of  high 
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torque  AMED’s  (such  as  reaction  wheels  and  CMG’s)  is  inconsistent  with  (or  at  least  tech¬ 
nically  very  challenging  for)  the  state-of-the-art.  Whereas  reaction  jets  can — depending  on 
physical  placement  issues — generate  very  large  torques  with  rapid,  discontinuous  action;  but 
their  use  in  throttled  (continuous)  mode  is  technically  challenging.  The  decomposition  of 
control  authority  suggested  in  the  following  design  approach  is  chosen  to  exploit  the  benefits 
of  current  actuation  technology  for  attitude  control. 

One  design  approach  for  rapid  slewing  focuses  on  controlling  the  dynamics  of  the  rigid 
body  motion.  We  discuss  the  implications  of  dynamic  reaction  from  the  elastic  appendage 
in  the  context  of  robust  control  performance.  Recall  that  the  rigid  body  dynamics  can  be 
given  in  the  form. 


7  =  r(7K  (4.49) 

d'i,  =  Ib  l(Ibu>b  x  u>b)  +  Ib  lrb  (4.50) 

where  y,u>b  G  Jf3  are  respectively  the  vector  attitude  parameters  and  angular  rates  referenced 
to  the  body  fixed  frame. 

To  control  the  vehicle  attitude  we  now  design  an  implementation  of  PLF  compensation 
for  the  system  outputs, 

V  =  7, 

given  the  external  torques,  rb,  applied  to  the  rigid  body.  We  remark  that  the  assumptions 
leading  to  the  model  (4.50)  together  with  the  choice  of  primary  system  outputs,  y,  indicate 
that  the  relative  degree  from  each  input /output,  pair,  rbi  >— *  y,  is  2.  Transformations  obtained 
as  in  (2.27)  indicate  that  the  system  can  be  exactly  linearized  by  feedback  (i.e.,  no  zero 
dynamics  are  obtained  in  the  transformation  to  the  normal  form). 

The  computation  of  the  normal  form  coordinates  is  obtained  as  in  section  2  by  differen¬ 
tiation  of  the  output.  In  this  case, 


-i  =  y  =  7  (4.51) 

z2  =  y  =  r(7)w  (4.52) 

and  since  the  3x3  matrix  T  is  nonsingular  for  almost  all  7  the  transformation  is  invertible. 
A  straightforward  computation  obtains  the  normal  form  as, 


~2  —  -4( z2)  +  B{z\,  z2)rb 

(4.53) 

where 

5(c1)  =  r(r1)/,r1 

(4.54) 

.4(^2)  =  [7r^r(7)w  +  n7 )ib\h^  x  u,)]| 

(4.55) 

Thus  the  feedback  compensation  which  linearizes  the  system  (4.50)  for  large 
is  [Dwy84], 

angle  motion 

Tb  =  fl(z,)-,{a(t)-^i,;:2)} 

=  -IbivujTy  -  Ibu)  x  +  IbT(y)a(t) 

(4.56) 
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where  a  6  9i3  is  the  “synthetic  control”  (in  acceleration  coordinates).  To  implement  the 
rapid,  large  angle  slewing  control  we  apply  the  bang-bang,  time-optimal  switching  logic  to 
each  axis;  i.e., 

—  ^inax  { /(  ~2,i)  “f-  2araaxZi  {} 

as  in  (4.32)  with  /(•)  as  in  (4.38).  This  form  assumes  identical  acceleration  limiting  is  desired 
in  each  axis. 

A  straightforward  computation  reveals  that  we  can  write, 


CN 

M 

>7 

ST 

ii 

n 

N 

>7 

(4.57) 

with 

d(-i,  z2)  =  [/yTu)Ii>u>  +  If,u>  x  w]  . 

,t'= *i .  u=r(ii  )i2 

(4.58) 

The  normal  form  equations  can  be  written, 

Z\  =  ~2 

(4.59) 

-2  —  -2  ){<^(  ^1 ,  Z2)  +  rb}- 

(4.60) 

Note  that  the  factorization  obtained  in  (4.57)  indicates  that  the  matching  condition  for 
robust  feedback  linearization  is  valid  for  the  model  terms  contributing  to  d(zi,z2). 

One  approach  to  implementation  of  feedback  linearizing  control  for  rapid  slewing  is  to 
implement  the  compensation  in  two  components; 


n  =  Tc  +  r\  (4.61) 

where  the  continuous  mode  torque  component  is, 

T  d{  ,  —2  )  I  -J  =r(7)u> 

=  -{■yT0jIl,uj  +  Ib<jJ  x  u;].  (4.62) 

This  component  compensates  by  decoupling  the  Coriolis  and  centripetal  accelerations  arising 
from  large  angle,  rapid,  rigid  body  rotations.  Note  a  typical  source  of  model  uncertainty 
arises  due  to  asymmetries  in  system  inertia  properties  which  are  typically  not  considered  in 
control  system  design  studies.  Such  asymmetries  arise  in  the  SBL  system  with  the  addition 
of  trackers  or  additional  mirror  components.  We  see  that  the  smooth  control  is  not  directly 
driven  by  the  subsequent  synthetic  control  which  obtains  the  desired  slewing  maneuver. 
Thus  the  smooth  control  functions  to  reliably  achieve  linearization. 

The  second  component,  r*,  will  be  implemented  by  discontinuous  or  saturation  mode 
actuation.  We  shall  now  indicate  a  simple  approach  to  identify  a  switching  control  law  for 
ts  based  on  VS  control  design.  After  application  of  the  control  (4.61)  to  (4.60)  we  obtain 
the  design  system  for  the  choice  of  the  saturation  mode  control  law, 


zi  —  ‘•i 

-  =  B(zi,  Z2)t* 


-2 


(4.63) 
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For  design  of  VS  control  we  follow  the  method  described  in  section  2.  Thus  we  choose 
the  switching  surface  in  the  normal  form,  ^-coordinates.  The  choice  of  the  switching  surface 
establishes  the  ideal  dynamic  behavior  in  sliding  and  its  design  should  embody  the  system 
design  goals.  In  section  2  and  [KK89]  we  outlined  options  for  stabilization  of  the  origin.  In  the 
present  design  the  large  angle  slewing  motion  is  critical  and  the  design  goal  is  time  optimal 
motion.  Solution  of  the  general  problem  of  time-optimal  attitude  control  for  the  rigid  body 
equations  (4.50)  with  independent  saturation  constraints  on  the  available  torques,  |r,|  <  rm(vx 
is  an  open  problem  [AF66].  The  formal  time-optimal  control  problem  does  not  embody 
considerations  for  robustness  or  implementation  issues  as  described  above  and  we  will  not 
consider  it  further.  Our  approach  is  to  develop  requirements  for  the  torque  authority  of  the 
saturation  mode  control,  r* ,  so  that  decoupled,  independent  single  axis  slewing  is  obtained. 

In  section  2  we  showed  that,  in  sliding,  the  inverse  dynamic  control  which  linearizes  the 
system  response  is  obtained  implicitly.  The  PLF  also  obtains  decoupling  of  the  individual 
response  of  the  motion  about  each  axis.  Our  approach  is  to  choose  the  switching  surfaces  for 
independent  axis  slewing  by  applying  the  modified,  quadratic  switching  surfaces  for  time- 
optimal  servocontrol  (4.31)  to  each  axis. 

To  make  these  statements  precise  we  introduce  the  following  notation.  In  the  normal 
form  coordinates  we  decompose  the  positions  and  rates  in  terms  of  decoupled  axis  motions 
as, 

zk,i  N 
~k,2 
~k,3  / 

for  k  =  1,2.  Then  in  the  normal  form  coordinates  we  establish  three  switching  surfaces, 

•*<(-)  =  f(z2,i)  4-  2amRX 
for  i  =  1,. . .  ,3.  Let  s(z)  -  [.Si,  s2,  53]T. 

It  is  easy  to  see  that  if  we  can  enforce  the  sliding  mode  over  a  significant  region  of  the 
switching  manifold,  M.,  then  the  independent  axis,  time-optimal  slewing  will  be  obtained. 
The  design  of  the  required  switching  control  to  achieve  sliding  on  M.,  and  to  guarantee 
reaching  the  manifold  from  the  expected  region  of  operation  after  a  sufficiently  small  transient 
can  be  based  on  Lyapunov  stability  arguments  as  follows. 

Consider  a  V-function, 

v  =  i»rUM=) 

defined  in  the  normal  form  coordinates.  The  time  derivative  of  V  along  trajectories  of  (4.63) 
is, 

V  =  sT(~-)i(~')  = 

=  sr(z)[2amRX/3,  F'(z)] 

=  *T(=)  {2amRXz2  +  F'(z)B(zl)r>}  ,  (4.64) 

(4.65) 


BW 


where 


F'{z2)  :=diag {/'(z2,.),  *  =  1,2,3} 
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and 

2^,  for  £  >  g2j2 

fit)  ■=  92,  for  |£|  <  g2/ 2  .  (4.66) 

.  — 2£,  for  (  <  -g2/2 

To  guarantee  that  the  switching  manifold  is  reached  in  finite  time  and  that  sliding  is 
achieved  we  seek  the  control  law  which  satisfies, 


min  sr(c)F,(;)B(z)rs  <  —  /.i 
r-e  T.„, 

for  f.L  >  0  with  admissible  controls 

Tad  =  {r  €  5?3  :  M  <  TmBX,  i  =  1,2,3}. 

Claim  4:  The  saturation  mode  control  law  which  is  optimal  for  ( 4.67 )— (4.68)  is, 


r,1  =  ~Tmnx  sgn{BT(z1){F'(=2)}Ts(z)} 


(4.67) 


(4.68) 


(4.69) 


for  i  =  1,2,3.  Together  with  the  optimal  switching  strategy  (4.69)  a  sufficient  condition  for 
the  required  saturation  mode  control  authority  to  guarantee  that  V  <  0  in  a  given  closed 
set,  7  €  B-,  C  9?3,  containing  the  origin  is 


Tmax  >  amBX  max,=li2,3  max^g^  1 5Z?=i [r(-)r)/6  1 


(4.70) 


Proof:  From  (4.64)  it  is  sufficient  to  show  that 


TmKX  >  max  max 


2a 


max  |  — '2  I 


11UVA  IlHtA  .  _  a  ,  _  _  .  ■ 

-1-2-3  l£j„,[F'(z2)B(*i)]y| 

From  (4.65)-(4.66)  and  (4.54)  we  can  simplify  the  ratio  for  each  i  as, 

2c*mBx|-2,i|  2a  max  |  *2  ,i  | 


I  E’=,[S'(Z2)S(..1)]„|  I  E‘=1  /'(z2,)(r(z,)/,-'|yf 
To  find  the  worst  case  for  each  i  from  the  discontinuous  definition  (4.66)  we  note  that  for 

I-2..I  >  W2  , 

2am«x  I  "2,t  I  2a  max  l^2,t| 


\f'M  I 


2|*,,.-| 


=  a„ 


and  for  the  case  \z2d  I  <  92/2  , 


2a 


max 


max  | ~2,t| 


2  a 


max  I  ~2,i  I 


1-2..  l<ffa/2  |/'(z2i;)| 


9  2 


=  Ocv 


*2., 


a 
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Figure  4.3:  Coordinated  Rapid  Slewing  Control  Implementation 

To  obtain  the  required  strict  negative  requirement  V  <  —p  similar  relations  can  be  de¬ 
veloped  which  involve  assumptions  on  the  a  priori  maximal  rates,  \z2\  <  7mnx.  Our  principal 
concern  is  in  implementing  the  bang-bang  control  with  sufficient  control  authority  TmfiX  such 
that  reaction  of  the  flexible  appendage  response  will  be  effectively  decoupled  from  the  rigid 
body  attitude  so  that  time-optimal  slewing  can  be  obtained.  Note  that  the  control  law 

(4.69)  can  be  implemented  independently  of  state  measurements  associated  with  the  flexible 
appendage  deformation.  The  overall  control  structure  is  illustrated  in  Figure  4.3.  From 
the  above  analysis  we  see  that  robustness  to  model  uncertainty  can  be  obtained  if  sufficient 
control  authority  in  slewing  torque  r*  (in  particular  Tmrix  is  large  enough)  is  available  so  that 

(4.70)  is  satisfied  for  all  model  variations  [GP82], 

4.3  Coordinating  Control  for  System  Slewing  and  Multibody  Alignment 

Throughout  the  modeling  and  control  development  for  multibody  slewing  we  have  empha¬ 
sized  the  importance  of  establishing  a  body  fixed  coordinate  system  attached  to  a  system 
primary  body.  The  dominant  system  nonlinearities  arise  from  the  introduction  of  nonlinear 
kinematics  and  dynamics  of  large  angle  rotations.  We  have  implicitly  assumed  that  struc¬ 
tural  deformation  contributing  to  multibody  misalignment  dynamics  are  small  and  therefore 
linear  when  referenced  to  the  body  fixed  frame.  Our  approach  for  coordinating  multibody 
alignment  and  rapid  slewing  will  be  based  on  these  assumptions.  We  now  show  that  appli¬ 
cation  of  PLF  compensation  for  linearizing  and  decoupling  the  natural  nonlinear  response 
of  rapid,  large  angle  slewing  motions  of  the  system  primary  body  results  in  decoupling  the 
multibody  alignment  problem  from  the  slewing  problem. 

Referring  the  general  multibody  dynamics  model  (4.1 )— (4.2),  consider  now  the  special 
case  of  (5.47) — (5.52).  Let  q  =  ['yT,uT}T,  p  =  [u/r,ur]r,  and  u  =  [7r,£T]r.  Then  the  SBL 
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model  for  multiaxis  slewing 


developed  in  section  5.1  can  be 


written  in  the  form, 


Mn(q) 

MlM) 


Mu{q) 

M-z2{q) 


“b  \  Bn(q,p)  B12{q,p ) 

«  /  Bf2(q,p)  B22(q,p) 


A'n  (q,p)  A'u  (q,p) 

(  7  )  - 

.  Ku{q,p)  A'22(?,p) 

\  u  ) 

Gb 

0 


Tb- 


(4.71) 


We  focus  attention  on  the  regulation  of  the  primary  system  body  frame  attitude  and  take 
as  primary  system  outputs, 

y  =  [C\,0]<?  =  C\  7, 

where  we  assume  that  y  and  rb  are  3- vectors  and  the  3x3  matrix  C is  full  rank.  We  further 
assume  that  the  regulated  outputs  (note  that  these  outputs  do  not  necessarily  coincide  with 
sensor  locations)  are  collocated  with  the  primary  body  torques; 


Under  these  assumptions  it  is  by  now  easy  to  see  that,  the  zero  dynamics  resulting  from  PLF 
compensation  are  given  by, 

M22u  4-  B22ii  q-  K22u  =  0.  (4.72) 

Note  that  (4.72)  is  linear  time  invariant  equation  of  the  free  response  of  the  multibody 
alignment  and  structural  vibration  with 


7  =  u>  =  0. 

Stability  of  the  zero  dynamics  is  a  prerequisite  for  stabilization  using  PLF,  but  for  practical 
implementations  the  degree  of  stability  may  be  an  important  factor.  For  the  SBL  rapid 
slewing  and  pointing  problem  this  decomposition  of  control  problems  is  critical  to  system 
performance. 

Decomposition  of  Control  Authority  and  Stabilization  of  the  Zero  Dynamics. 

Consider  the  general  model  of  (4.1 )— ( 4.2)  modified  to  include  additional  controls, 

q  =  f (q)p  (4.73) 

M(q)p  +  B(q,p)p  +  K(q,p)q  =  G^n  +  (4.74) 

primary  control  alignment  control 

A  direct  computation  based  on  (4.4)  indicates  that  the  required  PLF  compensation  for 
the  output  (4.3)  includes  feedforward  from  the  structural  controls; 

n  =  {CTM~'Gh)-lv  -  ( CtM-lGh)~xCtM-'[{MD  -  B)p  -  Kq  +  Gaf,l  (4.75) 

and  the  resulting  zero  dynamics  are 


Af22ii  +  B22ii  -f-  K22u  —  Gafa. 


(4.76) 
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For  application  to  rapid  slewing  of  the  multibody  system  the  primary  system  output  will 
be  taken  to  be  y  —  (or  in  terms  of  the  Gibbs  vector  y  =  7).  Referring  to  the  equations  of 
motion  for  the  generic  multibody  system  derived  in  Section  5.1  we  see  that  the  zero  dynamics 
(4.76)  are  merely  the  linear  dynamics  of  the  flexible  structure  with  cantilevered  boundary 
conditions.  The  structure  controls  effect  multibody  alignment  and  vibration  suppression 
relative  to  the  system  body  fixed  frame. 

The  design  approach  we  suggest  is  the  integration  of  precision  multibody  alignment  and 
vibration  suppression  based  on  robust  stabilization  of  (4.76)  together  with  the  PLF  compen¬ 
sation  and  slewing  design  described  above.  We  remark  that  the  design  approach  is  predicated 
on  the  fact  that  multibody  alignment,  for  precision  pointing  and  tracking  is  required  at.  the 
termination  of  the  slewing  maneuver.  Stabilization  of  (4.76)  by  the  introduction  of  local  feed¬ 
back  will  effect  transient  responses  during  the  slewing  maneuver  but.  the  principal  concern 
is  stabilization  at  the  end  of  the  slew. 

Remark:  The  decomposition  for  coordination  of  multibody  alignment  with  rapid  slewing 
can  be  viewed  as  a  generalization  of  implementations  commonly  used  for  vibration  suppres¬ 
sion  in  flexible  space  structures.  As  described  in  [Jos89]  implementation  of  feedback  control 
for  active  vibration  damping  commonly  uses  pairs  of  actuators  and  sensors  symmetrically 
located  so  as  to  eliminate  sensitivity  to  rigid  body  modes.  Such  configurations  are  required 
since  the  low  levels  of  control  authority  typically  available  for  structural  damping  are  too 
small  to  effect  attitude  and  translation  of  the  space  vehicle.  If  the  response  of  the  rigid 
body  modes  are  not.  decoupled  from  these  control  loops  the  actuators  would  saturate  during 
any  maneuver.  Although  the  approach  based  on  PLF  for  decoupling  clarifies  the  role  of 
local  linear  mode  vibration  damping  in  the  context,  of  the  nonlinear  system  dynamics.  Note 
however  that  in  general,  stabilization  of  the  local  linearization  of  the  zero  dynamics  may  not 
guarantee  global  stabilization  with  feedback  linearizing  control  [KS89] . 
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5  Simulation  Tradeoff  Studies  for  Multiaxis  Slewing  of  SBL  Sys¬ 
tem  Model 

In  this  section  we  develop  a  model  for  generic  SBL  system  beam  expander  as  a  two  body 
system  with  elastic  support  structure.  The  model  is  refined  from  that  developed  and  reported 
in  [BBKA88]  by: 

•  3  rotational  degrees-of-freedom  included  in  system  primary  rigid  body  motion, 

•  area  moments  of  attached  structure  modified  to  represent  tripod  metering  truss, 

•  system  inertias  and  slewing  requirements  scaled  to  represent  capabilities  of  ASTREX 
facility  and  initial  test  article. 

5.1  Multibody  System  Model  for  Multiaxis  Slewing  and  Precision  Alignment 

The  modeling  method  of  Lagrange's  equations,  as  described  in  the  report  [BBKA88],  along 
with  spatial  discretization  via  collocation  by  splines  to  develop  a  finite  dimensional  model 
suitable  for  simulation  of  large  angle,  multiaxis  motion  of  a  generic,  two-body  model  of  a 
SBL  system  beam  expander  is  described  in  this  section.  Beginning  with  the  generic  flexible 
spacecraft  model  [BBKA88,  sec  4]  the  following  model  assumptions  are  made: 

1.  axial  appendage  deformations  are  negligible,  r/3  %  r, 

2.  translation  velocity  of  the  system  primary  body  is  negligible,  R  0, 

3.  torsional  appendage  deformations  are  negligible,  V’(  - )  ^  0. 

Thus,  following  the  DPS  modeling  approach  discussed  in  [BBKA88,  sec  3],  the  configuration 
space  for  the  distributed  parameter  model  is  50(3)  x  Hg,  where  Hg  is  the  set  of  continuously 
differentiable  functions,  \(c)  =  [771,  tj2,  9,  d>}T  £  Hg  defined  on  the  interval  z  £  [0,£],  and 
which  satisfy  the  geometric  boundary  conditions  : 

Q  :  *71  (-)  =  0,  r U(z)  =  0.  9(z)  =  0,  <j>(z)  =  0,  at  z  =  0. 

With  the  notational  conventions  in  [BBKA88],  the  model  will  be  described  by  the  appendage 
lateral  deformation,  rj(t,z)  =  [771,772]  ,  along  the  body-fixed  x  and  y  axes,  respectively,  and 
appendage  angular  deformation,  £(c,£)  =  [9(z),4>{z)}T .  The  attitude  angular  rates  referenced 
to  the  body  fixed  frame  are  denoted,  u 

The  potential  functions  for  the  variational  analysis  leading  to  the  dynamic  equations  of 
motion  reduce  in  this  case  to: 


=  \u>l  IhUJh 

+  k  [  llfi67(~)  +  *7t(-)||2M  dz 
Jo 

+\ /oV*>  +  p(t(=)}Tn^  +  PM=)]p 


T 


(5.1) 
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Notation 

Explanation 

P 

mass  density 

A 

cross  section  area 

E 

elasticity 

kG 

I 

area  moment  of  inertia 

effective  shear  modulus 

•  dx 

*  =  Tt 

time  differentiation 

*=«•  =)=ff(M) 

partial  differentiation 

Table  5.1:  Standard  Notation  for  Lagrangian  Mechanics 


v(v,Z) 


i  ft{EIffe2z(=)  +  EI^l(z)  (5.2) 

Jo 

+  K!G>1[(77i(z));  —  #(~)]2  +  K2G.4[(r/2(z));  —  4>{z )]2}  dz,  (5.3) 


where 


mi)  =  J  J0  {nf-lVt  +  (f-2 (t  +  {r]tz)T -3Vtz  +  {Ztz)TZ4{Ztz)}dz, 


p  = 


■  0  0 
0  1 
1  0 


1  • 
0 
0 


(5.4) 


T  is  the  kinetic  energy,  V  the  potential  energy,  and  R  the  Rayleigh  dissipation  function.  The 
notation  used  is  standard  in  continuum  mechanics  of  beams  and  is  summarized  in  Table  5.1. 
The  area  moment  tensor  is  assumed  to  have  the  form2, 


/(;) 


/*«  0  0  ‘ 
0  Im  0 
0  0  J„  . 


For  simulation  studies  contained  in  this  report,  we  utilize  a  simple  damping  model  given 
by  the  assumption  that  the  matrix  coefficients  in  the  dissipation  function  are  of  the  form, 

—  _  Ci,«  0 

“  0  C2,  .  ' 

Note  that  Ei,E2  model  external  dissipation  while  Es.E^  model  internal  dissipation  [BK89J . 

The  control  forces  acting  on  the  system  primary  body  are  modeled  as  external  torques 
and  enter  the  variational  setup  through  a  virtual  work  expression  of  the  form, 

8W  = 


2 We  use  the  so-called  NASA  standard  or  321  convention  [Gol82] . 
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In  the  sequel,  we  refer  to  the  kinetic  energy  expression  in  terms  of  three  components; 


■^rigid  body 

=  Ib^h, 

(5.5) 

■^flex! 

=  \[  1  \^bV{:)  +  Vt(z)\\2pAdz, 

Jo 

(5.6) 

^flex2 

=  ?  Jo  [“b  +  tt(z)]T I[ub  +  it{z)]p  dz. 

(5.7) 

FEM:  Collocation  by  Splines 

A  model  based  on  the  Finite  Element  Method  (FEM)  is  obtained  by  spatial  discretization 
via  collocation  by  splines.  The  use  of  splines  for  such  purposes  is  described  in  [BBKA88]. 
Given  the  geometric  boundary  conditions  described  above  the  discrete  approximations  of  the 
scalar  valued  functions,  t/i ( ~ ) ,  #(-),  d>(~),  are  decoupled  and  reduce  to  approximation 

of  a  single  scalar  function,  say  7 (z),  on  the  interval  ;  6  [0,  f]  with  boundary  condition, 
7(0,  t)  =  0.  Dividing  the  interval  [0,  f]  into  N  uniform  subintervals  and  using  first-order 
B-splines  we  obtain  the  approximation: 

(5-8) 

*=o 

with  the  boundary  condition, 

N 

£  W)B'i £(0)  =  o. 

i= 0 

Then  following  the  reduction  procedure  described  in  [BBKA88,  §4.4.2]  we  obtain  the  FEM 
approximations  with  $(::)  as  given  in  [BBKA88,  Eqn.  (4.47)], 


Vi  (s,0  ~  $r(-)f?i(<) 

(5.9) 

V2  (-,*)  %  $T(z)ij2(t) 

(5.10) 

9{z,t)  %  $T{z)0{t) 

(5.11) 

%  $r(z)<£(f). 

(5.12) 

5.1.2 

Reduction  of  the  Kinetic  Energy  Function 

Using  the  FEM  approximations  of  order  N  we  obtain  the  spatial  discretization  of  the  kinetic 

energy  terms, 

^fleXi  ~  +  2Wj ,  4-  T]  Irjri^l} 

(5.13) 

7flex2  35  JutZ  +  Z 

(5.14) 

where 

'  vl +  <7  -ijf  Nnr}rj2 

Lu>{v)  =  2  ~vl  N, 7 

-tfNr,  +  vl _ 

(5.15) 
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N„  =  t  PA*(z)*T(z)dz, 
Jo 

<  =  j‘pAz*T{z)iz, 


Jujrj 


“'ll 


Ju, 


and 


— 


Jtt  = 


/ 

1  pAz2 

dz  = 

pAl 

Jo 

1 

3 

0 

-NT  1 

Nj 

0 

*?2  ^00 

if  Nnn  J 

0 

0 

7 

rl 

:  pl(z)dz 

7 

JO 

'  0 

^1 

- 

Nj 

0 

7 

0 

0 

'  Nee 

0  i 

0 

N<txj>  J 

The  following  terms  are  used  in  the  above  expressions, 


N, 


4> 

Nj 


Nj 


=  [l  ph;<t>T(z)dz, 

Jo 

=  [  />4y$r(z) 

Jo 

=  [‘ pIxx$T(z)dz, 

Jo 

h*  =  f  pIxx*{z)*T{z)dz, 

Jo 

Nee  =  f  pIyy$(z)<i>T(z)dz. 

Jo 

Note  that  the  total  kinetic  energy  can  be  written  in  the  standard  form, 


Ns 


/  U>t,  \ 

ij 

\  (  J 


where 


Mub,v,a  = 


lb  ~b  lutuiiv)  "b  Jw 

tljl) 


JL 


^00  0 


0  J, 


a 
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(5.16) 

(5.17) 

(5.18) 

(5.19) 

(5.20) 

(5.21) 

(5.22) 

(5.23) 

(5.24) 

(5.25) 

(5.26) 

(5.27) 

(5.28) 


(5.29) 


(5.30) 
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5.1.3  Reduction  of  the  Potential  Energy  and  Dissipation  Functions 

As  above,  the  potential  energy  and  dissipation  functions  are  reduced  by  substitution  of  the 
FEM  approximations.  For  the  potential  energy  we  obtain, 

V  =  ${9T  Ke9  +  <j>T  +  r]J  hmiji  +  fjT  hn2T}2  +  29T  +  24>T  (5.31) 


where 


Kg  = 

z)$r(s)  +  K  1GA$(z)<t>T(z)}dz 

(5.32) 

iu  =  I\eixx*z{ 

J  0 

z)^(-)  +  K2GA$(z)$T(z)}dz 

(5.33) 

A',,  =  j\^GAiMiTAz))iz 

(5.34) 

Km  =  [\k2GA$z 
Jo 

(*)«r  (-IK- 

(5.35) 

Kgm  =  I' {K-iGA${z)$*(z)}dz 

Jo 

(5.36) 

A’*a  =  -J  fl{K2GA*(z)*TA:)}dz. 

Jo 

(5.37) 

Similarly,  the  dissipation  function  reduces  to  the  form, 

R  —  |{r/!  Bmfh  +  jj2  Bl2fj2  +  9  Bg9  4-  4>  B^  +  29  Bgmfjl  +  2 4>  Bfafy}, 

(5.38) 

where  given  the  expressions, 

ft,  j 

f  <t>(z)*T(z)dz, 

0 

(5.39) 

K  -  j 

f‘*t{z)*T(z)dz, 

0 

(5.40) 

«  =  J 

f($z(z)<S>T(z)dz, 

0 

(5.41) 

we  find  it  convenient  to  express  the  individual  matrix  coefficients  in  R  (arising 

from  the 

FEM  approximation)  as: 

Br,y  -  Cil  El 

,  +  C>3  R'o  for  1  =  1,2, 

(5.42) 

B$  =  (Cl2  +  Cl3)-Ro  +  CuR'0, 

(5.43) 

B<t,  =  (C22  +  (2z)Ro  +  (24R0, 

(5.44) 

Bme  =  -Ci3#o. 

(5.45) 

Bm<t>  =  -£23  R'o- 

(5.46) 

Finally,  we  summarize  the  model  stiffness  and  damping  using  the  expressions, 


>, 

•4 

II 

Km 

0 

0 

K  rj2 

43 

^>v 

II 

R  9r) j 

0 

0 

A  tpTJj 

>, 

/>* 

II 

'  Kg 
0 

0 

K* 

Bn  = 

‘  Bm 
0 

0 

Bfn  . 

1  Bnt  — 

B«n  1 

0 

0 

B<t>  m 

.  *<  = 

Bg 

0 

O  qq* 
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5.1.4  Lagrange’s  Equations 

Finally,  we  compute  the  equations  of  motion  using  Lagrange’s  equations  in  the  form: 


<?  = 
M(q)p  + 


f(<7)p 


dM(q)p 

1 

\dM(q)p] 

.  d<l  . 

P  -  2 

dq  . 

SV  6R  „ 
p  +  Tq-~Jp^Q' 


(5.47) 

(5.48) 


with  generalized  coordinates  q  —  {(,b,V,Z},  and  velocities,  p  =  { ,  7? ,  ^ } . 


Equations  of  Motion  for  Multiaxis  Slewing  Applying  Lagrange’s  equation  results 
(after  some  simplification)  in  the  equations  of  motion  with  kinematics  of  the  system  attitude 
expressed  in  terms  of  the  Gibbs  vector,  7; 

7  =  r(7)ufc,  (5.49) 


and  kinetics, 


[/&  "I"  Iw  "b  J u ~b  luiriP  +  "b  “b  luJW  “b  "b  luriq  ~b  “b  H TJ  —  Tb  (5.50) 

+  KnV  ~  \HT ^b  +  Rfj  4-  Knfj  +  Knit  +  Bnq  +  B^l  =  Gnfs  (5.51) 

J&b  +  jJ  +  Kti  +  K%f}  +  B<£  +  =  G{T,  (5.52) 


where 

H(fj,ub,v) 


-“>2  il  N*,  -  U/3N* 

2 V2 

-U1V2  Nnv  +  2u}^\  Nm 

-uivl  Nnv  ~  N* 

.  (5.53) 

-Wi  N*  -b  u>3^vl Nm  +  Nvn 

+  <M32fj2  Nnr]  +  fj 2  Nm  . 

RM  =  °w 

0 

(5.54) 

5.2  Computer  Simulation  Model  for  Nonlinear  Control  Law  Tradeoff  Studies 

A  computer  simulation  model  was  developed  and  coded  based  on  the  above  FEM  model  for 
attitude  slewing  of  a  flexible  space  craft  ( 5.49)— ( 5.52 ) .  The  model  parameters  and  geometry 
was  chosen  to  represent  the  dynamics  of  a  SBL  system  beam  expander  with  structural 
flexure  arising  from  the  flexure  response  of  the  metering  truss.  A  goal  of  the  second  year 
effort  was  to  develop  concepts  for  experimental  testing  of  nonlinear  control  concepts  for 
flexible  space  structures.  A  suitable  facility  for  testing  large  angle,  multiaxis  slewing  and 
precision  pointing  is  currently  under  construction  at  the  Astronautics  Laboratory.  The 
ASTREX  facility  with  include  3  axis  slewing  and  the  initial  test  article  is  a  scaled  model 
of  a  SBL  beam  expander.  We  have  adapted  our  nonlinear  slewing  control  simulation  model 
to  represent  available  parameters  and  geometry  of  the  ASTREX  test  article.  Details  of  the 
modeling  assumptions  and  geometry  are  given  below. 
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5.2.1  Simulation  model  geometry  assumptions  and  parameters 

The  objective  of  the  simulation  model  design  was  to  generate  a  simplified  model  that  would 
duplicate,  as  closely  as  possible,  the  geometry  and  dynamics  of  the  ASTREX  test  article. 
Information  describing  the  test  article  was  sparse,  therefore,  the  following  assumptions  were 
required  to  develop  a  working  simulation  model: 

1.  The  basic  configuration  of  the  model  consisted  of  a  rigid  body,  representing  the  main 
mirror  assembly  and  associated  structure,  and  an  attached  tripod  structure,  represent¬ 
ing  the  flexible  metering  truss  with  secondary  mirror,  see  Figure  5.1. 


MET  =  2!M(C 
TRUSS 


Figure  5.1:  Simulation  model. 


2.  The  metering  truss  geometry  consisted  of  a  tripod  of  tubes  located  at  the  corners  of 
an  equilateral  triangular  base.  The  generic  model  discussed  above  was  readily  adapted 
to  represent  this  geometry  by  including  spatial  variation  of  the  area  moments  Ixx(z ) 
and  Iyy(z)  based  on  this  geometry.  Details  of  the  computation  of  the  area  moments 
and  the  resulting  FEM  model  approximations  are  given  in  Appendix  A. 


3.  The  dimensions  of  the  metering  truss  model  were  as  follows: 


Length, 1 
Base  width,  b 
Tube  inner  dia. 

Tube  outer  dia. 

Mass  of  secondary  mirror 


5.0  m 
3.8  m 
8.0  cm 
12.7  cm 

190  kg  (5  %  of  total  sys.  mass) 


4.  The  overall  system  characteristics  were  as  follows: 
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Mass  moments  of  inertia 

/» 

=  16640  kg  m2 

lyy 

=  16590  kg  m2 

/-. 

=  11660  kg  m2 

System  mass 

m 

=  3810  kg 

First  modal  freq.  of  metering  truss 

=  12.2  Hz 

5.  The  dynamics  of  the  secondary  mirror  were  neglected  and  it  was  modeled  as  a  point 
mass  at  the  end  of  the  truss. 

6.  The  modulus  of  elasticity  of  the  truss  material  was  chosen  so  that  its  first,  modal 
frequency  matched  that  of  the  ASTREX  test  article. 

5.3  Simulation  Results  for  PLF  Control  with  Rapid  Slewing  and  Precision 
Pointing 

Based  on  the  above  simulation  model  a  digital  simulation  was  developed  using  PC'-Matlab 
to  perform  tradeoff  studies  for  implementation  rapid  slewing  and  precision  pointing  of  an 
SBL  beam  expander  with  PLF  compensation. 

The  model  equations  ( 5.49)— ( 5.52 )  with  position  coordinates,  q  —  [7^,  t/T,  £r]r,  and  ve- 

T  -JT 

locities,  p  =  ]T<  can  be  reduced  to  the  generic  form  (4.1 )— (4.2).  The  simulation 

models  and  parameters  assumed  for  the  tradeoff  studies  are  summarized  in  Table  5.2.  Sim¬ 
ulated  time  responses  for  each  model  with  several  different  control  laws  are  displayed  in 
the  figures  contained  in  Appendix  C.  Each  simulated  maneuver  is  summarized  in  a  set  of 
figures  including  the  responses  of  position  7,  attitude  rates,  7,  system  torques  applied  to  the 
primary  body,  rh,  and  deformations  at  the  apex  of  the  truss,  77(f).  £(f).  Table  5.3  describes 
the  simulation  runs  performed  and  illustrates  the  relationship  between  system  models  and 
control  laws  used. 

Slewing  with  Smooth  PLF  Compensation.  For  precision  pointing  of  the  system  at¬ 
titude  we  take  the  primary  system  outputs  as  y  =  •)},.  PLF  compensation  for  the  resulting 
system  is  given  by  the  formula,  (4.4).  Multiaxis  slewing  and  pointing  control  was  imple¬ 
mented  on  an  independent  axis  basis  subsequent  to  PLF  compensation.  Thus  the  PLF 
Slewing  control  is  obtained  by  setting  the  v  controls  appearing  in  (4.4)  as, 

r,  =  cti,  for  t  =  1,  2, 3 

where  the  commanded  accelerations,  a;  are  generated  by  the  independent  axis,  time-optimal 
slewing  control, 

(at(f)  —  —  <a,nox  sat{/(7; )  -)-  2omRX7i }  1 

as  described  in  Section  4.1. 

Simulation  tradeoff  studies  completed  this  year  focused  on  the  implementation  of  PLF 
compensation  for  various  reduced-order  models  of  the  metering  truss  flexure.  The  reduction 
methods  described  in  Section  3  were  used  to  obtain  implementation  alternatives  for  PLF 
compensation  in  terms  of  the  quasi-rigid  model  with  DOF,  r  =  3,4,5,11.  The  case  r  =  11 
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Parameter 


Simulation  Model 


MODI  MOD2  MOD3  MOD4  MOD5  MOD6  MOD"  MOD8  MOD9 


p  -  mass  density 
E  -  elasticity 
N  -  #  finite  elements 
e  damping  factor 
TTTy  J  VTlfttf 


ODE  tolerance 


1520 

SBL  beam  expander  model  parameters 

200 

o 

- 

- 

- 

- 

- 

- 

- 

it 

.001 

.001 

.004 

.004 

.001  .001 

.001 

.001 

.001 

0.05 

0.0 

0.05 

0.0 

0.0  0.0 

0.05 

0.0 

0.05 

.349 

_ 

Control/ Actuator  parameters 

_ 

N/A 

N/A 

N/A 

N/A 

9.3  1.4 

27.0 

22. 

22. 

6.477 

6.477 

6.477 

6.477 

5.E6  5.E6 

5.E6 

5.E6 

5.E6 

0.0374 

- 

- 

- 

- 

- 

- 

- 

0.001 

- 

- 

Simulation  Solution  Precision 

- 

- 

Table  5.2:  Matrix  of  Simulation  Models  Considered 


Models 

Control  Laws 

Continuous  PLF/Slewing 

Disontinuous  PLF/Slewing 

a 

a 

a 

a 

Slew  rate 
limited  torque 

Discontinuous 

torque 

MOD4 

MOD5 

MOD6 

MOD7 

MOD8 

MOD9 


Multiaxis  Slew 


SIM11  SIM12  SIM  13  S 
SIM21  SIM22  SIM23  S 
SIM31  SIM32  SIM33 
SIM41  SIM42  SIM43 


SIM  5 
SIM6 


SIM  15 
SIM25 


SIM8 

SIM9 


Multiaxis  Large  Angle  Slew 


SIM16 

SIM26 


Table  5.3:  Matrix  of  PLF/Slewing  Control  Implementations 
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Slew  Times  (sec)  to  I7J  <  .001 


PLF  Model  DOF 

T  = 

3 

4 

5 

11 

Model  MODI 

Axis  1 

.8005 

.8218 

.8268 

.8258 

2 

.6709 

.6725 

.6933 

.6930 

3 

.5145 

.5162 

.5241 

.5223 

Model  MOD2 

Axis  1 

.8257 

.8246 

.8296 

2 

Mm 

.8335 

.6868 

.6872 

3 

ns 

.5168 

.5189 

.5183 

Model  i 

V10D3 

Axis  1 

>  1.00 

.8260 

.8273 

2 

.8497 

.8494 

.6892 

3 

.5172 

.5193 

.5227 

Model  ] 

V10D4 

Axis  1 

.8094 

.8277 

.8302 

2 

.6732 

.6920 

.6888 

3 

.5184 

.5136 

.5191 

Table  5.4:  Slewing  Times  for  Standard  3-axis  Maneuver 

corresponds  to  compensation  for  all  the  simulated  elastic  degrees  of  freedom  of  the  flexible 
structure  response  and  was  included  as  a  baseline  for  ideal  PLF  compensation. 

Simulation  runs  SIM  11 — SIM43  were  made  to  compare  slewing  time  and  peak  torque 
requirements  for  implementation  of  PLF  compensation  for  multiaxis  slewing.  For  com¬ 
parison  of  subsequent  time  responses,  a  nominally  small  angle,  multiaxis  maneuver  was 
considered  with  desired  system  attitude  given  by  7  =  0  and  initial  orientation  given  by 
7  =  [.07,  .05,  .03]t.  This  attitude  maneuver  can  be  expressed  in  Euler  angles  as, 

[V\M  =  [1.4°, -1.7°, 2.3°]. 

Performance  results  for  continuous  mode  PLF  implementation  are  summarized  in  Ta¬ 
ble  5.4  and  5.5.  The  performance  sensitivity  to  the  presence  of  a  significant  mass  (approx¬ 
imately  5  %  of  the  system  mass)  at  the  apex  of  the  truss  was  also  examined.  Despite  the 
relative  stiffness  of  the  structure  the  multibody  coupling  is  evident  in  the  attitude  overshoot 
r^r  the  small  angle  maneuver.  (All  simulation  runs  are  included  in  Appendix  C.)  Imple¬ 
mentation  of  the  PLF  compensation  for  multiaxis  slewing  requires  increased  bandwidth  and 
peak  torques  on  the  system  rigid  body.  The  importance  of  damping  in  the  appendage  can  be 
seen  in  Tables  5. 4-5. 5.  With  increased  damping  peak  torque  required  for  PLF  compensation 
will  decrease. 

To  illustrate  the  large  angle  response  with  PLF  compensation  we  also  simulated  a  nia- 
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Peak  Torque  ATpk(r )/T'pfe( 3 ) 
Relative  to  rigid  body  case 

simulation 
model  r  = 

PLF  Model  DOF 

3 

4 

5 

11 

MODI 

0 

1.11 

1.10 

1.12 

M0D2 

0 

.278 

.282 

.293 

M0D3 

0 

.848 

.842 

M0D4 

0 

.238 

.243 

Table  5.5:  Relative  Peak  Torque  Increase  with  PLF  Compensation  DOF 

neuver  characterized  by  7  =  [.5,  .4,  .3]r,  which  can  be  expressed  in  Euler  angles  as, 

[V\M  =  [45.°, -19.5°,  61.3°]. 

The  simulated  time  responses  with  PLF  compensation  for  the  rigid  body  modes  only  reveals 
significant  overshoot  and  in  the  case  of  a  mass  at  the  apex  significant  torque  switching  activity 
is  generated.  The  addition  of  PLF  compensation  for  the  rigid  body  plus  first  two  elastic 
modes  of  the  beam  recovers  ideal  time-optimal  attitude  orientation.  This  also  noticeably 
reduces  chattering  activity  near  the  time-optimal  switching  surface.  For  these  large  angle 
maneuvers  (probably  difficult  to  test  experimentally  at  the  ASTREX  facility)  the  lateral 
deformation  of  the  apex  (with  tip  mass)  is  less  than  1  cm.  indicating  the  stiffness  of  this 
structure. 

Slewing  with  Discontinuous/Continuous  Actuation.  Simulation  studies  SIM51-SIM91 
were  performed  using  coordination  of  discontinuous  and  continuous  modes  of  actuation  for 
PLF  compensation  and  rapid  slewing  as  described  in  section  4.2.  The  simulation  trade¬ 
off  studies  for  this  case  focused  on  the  requirements  for  increased  peak  torque  to  recover 
near  ideal  time-optimal,  independent  axis  slewing  using  only  the  rigid  body  attitude  output 
information  for  control  law  implementation. 

As  discussed  in  Section  4.1  actuators  are  often  subject  to  slew  rate  limits  which  will 
limit  switching  response  for  on-line  implementation  of  discontinuous  mode  control  laws  of 
the  form  (4.69).  Simulation  runs  were  completed  with  slew  rate  limited  implementation  of 
the  discontinuous  controls,  r*,  by  replacing  the  sgn  function  with  a  saturation  function  with 
linear  slope,  g\  (cf.  Table  5.2). 

The  significance  of  implicit  PLF  compensation  achieved  in  sliding  is  apparent  by  com¬ 
parison  of  the  time  responses  of  SIM13  and  SIM72.  In  both  cases  the  models  include  a 
secondary  body  at  the  apex  of  the  metering  truss  which  involves  significant  multibody  dy¬ 
namic  response.  Using  continuous  mode  actuation,  PLF  compensation  is  maintained  during 
the  entire  slewing  transient  by  control  torques  applied  to  the  primary  body.  Using  discon¬ 
tinuous  actuation  PLF  compensation  is  achieved  after  an  initial  transient — once  the  sliding 
condition  is  attained.  Thus  during  the  acceleration  phase  of  system  slewing  the  action  of  the 
continuous  compensation  component,  rc  does  not  compensate  for  multibody  dynamics.  It 
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is  also  clear  that  increased  torque  authority,  Tmax,  is  necessary  to  maintain  sliding  over  that 
required  for  purely  continuous  mode  PLF  compensation  of  SIM13.  The  primary  difference  is 
that  the  continuous  mode  PLF  compensation  achieves  decoupling  by  ‘in-phase’  compensa¬ 
tion  of  the  multibody  motion  while  the  discontinuous/continuous  mode  approach  of  SIM72 
is  achieved  without  explicit  feedback  of  the  motion  of  the  secondary  body. 

It  is  not  clear  at  this  point  in  our  study  how  significant  it  is  to  maintain  PLF  compensation 
during  initial  transients.  For  system  slewing  the  total  slew  time  and  terminal  pointing 
accuracy  are  clearly  dominant  concerns.  However,  it  is  possible  to  include  compensation  for 
multibody  motion  in  the  combined  implementation  and  we  hope  to  investigate  this  further 
in  the  next  year. 

Our  simulation  studies  of  the  method  of  PLF  compensation  and  slewing  by  coordinated 
discontinuous  and  continuous  modes  of  actuation  should  be  viewed  as  preliminary.  A  goal 
of  the  next  years  activities  will  be  to  further  investigate  implementation  alternatives  and 
performance  tradeoffs  in  this  area.  The  simulation  results  obtained  to  date  reveal  that 
implicit  PLF  compensation  and  near  time  optimal  slewing  can  be  obtained  with  increased 
system  torques.  We  believe  that  the  peak  torque  requirements  for  slewing  can  be  relaxed 
using  dual  mode  actuation  if  the  secondary  body  motion  can  be  taken  into  account  in  the 
continuous  mode  component  rc  of  the  torque  acting  on  the  rigid  body. 
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6  Conclusions  and  Directions 

The  second  year  study  has  provided  some  new  insights  in  implementation  alternatives  for 
coordinated  multiple  actuation  mechanisms  for  achieving  rapid,  large  angle,  multiaxis  slewing 
control  of  multibody  systems.  The  simulation  studies  illustrate  several  important  design 
issues  in  implementation  of  PLF  control  for  slewing  and  pointing.  Implementation  of  PLF 
control  by  continuous  mode  actuation  enforces  decoupling  during  system  transient  responses 
by  additional  control  authority.  Its  implementation  will  require  available  measurements  of 
structural  flexure.  Implementation  using  discontinuous  mode  actuation,  based  on  sliding 
mode  control,  enforces  feedback  linearization  implicitly  and  only  after  an  initial  transient. 
An  open  question  is  the  importance  of  enforcing  feedback  linearization  during  such  transients 
for  retargeting  and  pointing  maneuvers. 

Simulation  studies  have  been  scaled  as  much  as  possible  based  on  available  data  and 
system  parameters  for  the  ASTREX  test  article.  It  is  hoped  that  the  control  concepts  for 
nonlinear  PLF  compensation  and  slewing  can  be  tested  experimentally  on  the  ASTREX  or 
similar  facility.  A  goal  of  the  third  year  effort  is  to  develop  designs  for  control  law  vali¬ 
dation  using  laboratory  simulations  and  experiments  for  the  control  concepts  developed  in 
this  study.  Although  detailed  designs  and  laboratory  simulations  are  outside  the  scope  of 
the  present  study  we  hope  to  describe  a  framework  for  an  experimental  test  plan  which  will 
identify  significant  nonlinear  dynamic  responses  and  validate  the  importance  of  feedback 
linearization  for  their  compensation.  We  feel  that  the  presence  of  such  nonlinear  dynamical 
interactions  can  be  an  important  limiting  factor  in  predicting  the  responses  of  control  system 
designs  for  multibody  systems.  It  is  widely  recognized  that  such  nonlinear  dynamics  will 
be  evident  in  rapid,  large  angle,  multiaxis  slewing.  We  feel  that  in  the  presence  of  signifi¬ 
cant  multibody  coupling  that  pointing  and  tracking  precision  may  also  be  limited  by  weak 
nonlinear  interactions.  If  such  limitations  are  significant  they  can  best  be  demonstrated  in 
laboratory  simulations. 


REFERENCES 

References 


62 


[AB87]  0.  Akhrif  and  G.L.  Blankenship.  Using  computer  algebra  in  the  design  of  non¬ 

linear  control  systems.  Proc.  ACC,  1987. 

[AB88]  0.  Akhrif  and  G.  L.  Blankenship.  Robust  stabilization  of  feedback  linearizable 

systems.  In  Proc.  27th  IEEE  C'ont.  Dec.  Conf.,  1988. 

[ABB89]  0.  Akhrif.  G.L.  Blankenship,  and  W.H.  Bennett.  Robust  control  for  rapid  reori¬ 

entation  of  flexible  structures.  In  Proc.  1989  ACC ,  pages  1142-1147,  1989. 

[AF66]  M.  Athans  and  P.  L.  Falb.  Optimal  Control.  McGraw-Hill,  1966. 

[Akh89]  0.  Akhrif.  Nonlinear  Adaptive  Control  with  Application  to  Flexible  Structures. 

PhD  thesis,  University  of  Maryland,  College  Park,  MD,  1989. 

[Arn78]  V.I.  Arnold.  Mathematical  Methods  of  Classical  Mechanics.  Springer- Verlag, 
New  York,  1978. 

[BBKA88]  W.  H.  Bennett,  G.  L.  Blankenship,  H.  G.  Kwatny,  and  0.  Akhrif.  Nonlinear 
dynamics  and  control  of  flexible  structures.  Technical  Report  SEI-88-11-15-WB, 
SEI,  November  1988.  AFOSR/AFSC  contract  F49620-87-C-0103. 

[BI84]  C.I.  Byrnes  and  A.  Isidori.  A  frequency  domain  philosophy  for  nonlinear  systems, 
with  application  to  stabilization  and  adaptive  control.  In  Proc.  IEEE  CDC ,  Las 
Vegas,  NV,  1984. 

[BI85]  C.  I.  Byrnes  and  A.  Isidori.  Global  feedback  stabilization  of  nonlinear  systems. 
In  Proc.  2Jfth  IEEE  CDC ,  pages  1031-1037,  Dec.  1985. 

[BK89]  W.H.  Bennett  and  H.G.  Kwatny.  Continuum  modeling  of  flexible  structures  with 
application  to  vibration  control.  AIAA  J.,  27(9):  1264—1273 ,  September  1989. 

[BL87]  J.  Baillieul  and  M.  Levi.  Rotational  elastic  dynamics.  Physica  27D,  27:43-62, 

1987. 

[Bro78]  R.W.  Brockett.  Feedback  invariants  for  nonlinear  systems.  In  Proc.  6th  IFAC 
World  Congress ,  pages  1115-1120,  Helsinki,  1978. 

[Dwy84]  T.  A.  W.  Dwyer.  Exact  nonlinear  control  C  rge  angle  rotational  maneuvers. 
IEEE  Trans,  on  Auto.  Control,  AC-29(9):76  1,  1984. 

[Dwy88]  T.A.W.  Dwyer.  Slew-induced  deformation  shaping.  In  Proc.  27th  IEEE  CDC, 

1988. 

[DZM88]  R.  A.  DeC'arlo,  S.  H.  Zak,  and  G.P.  Matthews.  Variable  structure  control  of 
nonlinear  multivariable  systems:  A  tutorial.  Proc.  IEEE,  76( 3 ) :2 1 2— 232 ,  1988. 


REFERENCES 


[FH87] 

[Fre75] 

[GL76] 

[Gol82] 

[GP82] 

[Hir79] 

[HSM83] 

[Isi85] 

[Jos89] 

[KC87J 

[KK89] 

[KK086] 

[K088] 

[Kra87] 

[Kre73] 

[KS89j 

[Le87] 


63 

R.  B.  Fernandez  and  J.  K.  Hedrick.  Control  of  multivariable  nonlinear  systems 
by  the  sliding  mode  method.  Int.  J.  Control ,  46(3):  1019—1040,  1987. 

E.  Freund.  The  structure  of  decoupled  nonlinear  systems.  Int.  J.  Control , 
21(3)443-450,  1975. 

S.  Gutman  and  G.  Leitman.  Stabilizing  control  for  linear  systems  with  bounded 
parameter  and  input  uncertainty.  In  Proc.  7th  IFIP  C'onf.  on  Optimization  Tech¬ 
niques ,  page  729.  Springer- Verlag,  1976. 

H.  Goldstein.  Classical  Mechanics.  Addison-Wesley,  Reading,  MA,  1982. 

S.  Gutman  and  Z.  Palmor.  Properties  of  min-max  controllers  in  uncertain  dy¬ 
namical  systems.  SIAM  J.  Control .  20(6 ):850— 861 ,  1982. 

R.  M.  Hirschorn.  Invert ibility  of  nonlinear  control  systems.  SIAM  J.  Optim.  and 
Control.  17(2  ):289— 297,  1979. 

R.  Hunt,  R.  Su,  and  G.  Meyer.  Global  transformation  of  nonlinear  systems. 
IEEE  Trans.  Automatic  Control.  AC-28:24-31,  1983. 

A.  Isidori.  Nonlinear  Control  Systems:  An  Introduction.  Springer- Verlag,  1985. 

S.  M.  Joshi.  Control  of  Large  Flexible  Space  Structures.  Springer- Verlag,  1989. 

C.  Kravaris  and  C'-B  Chung.  Nonlinear  state  feedback  synthesis  by  global  in¬ 
put/output  linearization.  AIC'hE  Journal ,  33( 4 ) : 592—603 .  April  1987. 

H.G.  Kwatny  and  H.  Kim.  Variable  structure  control  of  partially  linearizable 
dynamics.  In  Proc.  1989  ACC,  pages  1148-11 53,  1989. 

P.V.  Kokotovic,  H.  Khalil,  and  J.  O'Reilly.  Similar  Perturbation  Methods  in 
Control:  Analysis  and  Design.  Academic  Press,  New  York,  1986. 

F.  Khorrami  and  U.  Ozguner.  Singular  perturbation  analysis  of  a  distributed 
parameter  model  of  a  flexible  manipulators.  In  Proc.  1988  ACC ,  pages  1701- 
1709,  1988. 

C.  Kravaris.  On  the  internal  stability  and  robust  design  of  globally  linearizing 
control  systems.  In  Proc.  ACC ,  pages  270-279,  1987. 

A.  J.  Krener.  On  the  equivalence  of  control  systems  and  linearization  of  nonlinear 
systems.  SIAM  J.  C’ont.  and  Optim.,  11:670.  1973. 

P.  V.  Kokotovic  and  H.  J.  Sussinan.  A  postive  real  condition  for  global  stabiliza¬ 
tion  of  nonlinear  systems.  Syst.  C'ntrl.  Letters ,  13:1 25—  1 33 ,  1989. 

M.  Lock  and  et  al.  Structural  dynamic  response  of  a  space  based  laser  system. 
Technical  Report  AD-B109  212,  February  1987. 


TSI-89-12-12-WB 


64 


[MC80] 

[RR87] 

[SKK87] 

[SS83] 

[SV87] 

[Wer78] 


G.  Meyer  and  L.  C'icolani.  Applications  of  nonlinear  system  inverses  to  automatic 
flight  control  design-systems  concepts  and  flight  evaluations.  In  P.  Kent,  editor, 
Theory  and  Applications  of  Optimal  Control  in  Aerospace  Systems ,  pages  10—1— 
29.  NATO  AGARD,  1980.  AGARDograph  251. 

L.  C.  Rogers  and  K.  E.  Richards.  Pacoss  program  status  and  results.  Technical 
report,  Air  Force  Wright  Aeronautical  Laboratory,  1987. 

M.  Wr.  Spong,  K.  Khorasani,  and  P.  V.  Kokotovic.  An  integral  manifold  approach 
to  the  feedback  control  of  flexible  joint  robots.  IEEE  J.  Robotics  and  Auto .,  RA- 
3(4  ):29 1—299,  1987. 

J.J.  Slotine  and  S.S.  Sastry.  Tracking  control  of  nonlinear  systems  using  sliding 
surfaces  with  application  to  robot  manipulators.  Int.  J.  Control,  38(2):465— 492, 
1983. 

M.  W.  Spong  and  M.  Vidyasagar.  Robust  linear  compensator  design  for  nonlinear 
robotic  control.  IEEE  J.  Robotics  Auto.,  R A- 3( 4 ) : 345—3-5 1 ,  1987. 

J.R.  Wertz,  editor.  Spacecraft  Attitude  Determination  and  Control.  D.  Reidel, 
Dordrecht,  Holland,  1978. 


TSI-89-12-12-WB 


65 


A  Advantages  of  Gibbs  vector  description  of  nonlinear  kinemat¬ 
ics. 


Dwyer  [Dwy84]  uses  the  Gibbs  vector  parametrization  of  the  rigid  body  attitude  to  com¬ 
pute  an  exact  feedback  linearizing  transformation  for  rigid  body  attitude  dynamics.  The 
simple,  algebraic  form  of  the  resulting  kinematic  relations  aid  in  the  computation  of  PLF 
compensation.  In  this  appendix  we  review  the  basis  for  certain  relations  useful  in  the  general 
case. 

Euler’s  theorem  says  that  a  general  rotation  between  three  dimensional  orthogonal  frames 
can  be  described  by  rotation  through  an  angle  9  about  a  unit  vector  e.  Thus  any  L  £  50(3) 
can  be  given  by, 

L(e,8 )  =  cos  9I3  +  sin  9  f2(e)  +  ( 1  —  cos  9)eer 


where 


0 


Q(x)  :  — 


x3 


-x2 


~X3  X  2 

0  —  xi 

x\  0 


In  terms  of  the  Gibbs  vector  [YVer78,  pp.  416],  7  £  3?3,  given  by, 


7  :=  tan 


e, 


(A.l) 


(A.2) 


the  general  rotation  (or  direction  cosine  matrix  [Wer78])  can  be  expressed, 

T,  N  (1  -  lTl)h  +  277r  -  D(7) 

M7)  =  - 7~. — f - • 

1  +  7 3  7 

In  terms  of  the  Gibbs  vector  the  kinematic  relation  for  rigid  body  rotation  [Wer78,  pp. 
513]  can  be  expressed  as, 


7  =  r(7)u/ 

=  -[w  -  w  x  7  +  (u;  •  7)7], 

where, 

f  (7 )  =  \[h  +  77T  +  0(7)]- 


(A.3) 
(A. 4) 


Facts  about  skew  symmetric  matrices.  The  following  relations  summarize  the  salient 
geometric  properties  of  50(3).  Let  fi  :  J?3  — *  50(3)  (as  defined  above)  and  take  fi  as  the 
inverse  map.  Let  .4  £  50(3)  and  u,  v  £  5?3.  Then  the  following  relations  hold. 

1.  v  x  u  —  fi([D(v),  0(tt)])  3 

2.  v  x  u  =  (7(uvT  —  vuT) 

3.  Au  =  D(/t)  x  u 

4-  £(u  xv)  =  Q(-v) 


3Hete  [/l,  B]  -  AB  -  BA,  the  matrix  commutator. 
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Proof:  See  Baillieul  and  Levi  [BL87]  for  proofs  of  (i)-(iii).  To  show  (4)  note  that  from  (3) 


0  —  «3  111 

'  V\  \ 

(  U2V3  -  uzv2  \ 

U  X  V  = 

u3  0  -«1 

v2 

= 

W317  -  uxvz 

— u2  U\  0 

V  ^3  ) 

\  uxv2  -  U2v  1  / 

,  and  we  have, 


—~{u  x  v) 
du 


Q 

—  {Q(u)u} 

0  t»3  -t'2 

-t’3  0  Vi 

1’2  —  Vi  0 


(A. 5) 
(A. 6) 


□ 


Lemma:  Using  the  Gibbs  vector  (A.2)  the  following  kinematic  relations  hold; 


d 

—  (r(7)^}  =  \[~fTujI3  4-  7^r  +  fi(-w)], 


d-f 


{r(7)w}  r(7)w  =  7Tu>r(7)u;, 


Proof:  To  show  (A. 7)  compute  the  indicated  differentiation  of  (A. 3)  using  fact 
show  (A. 8)  use  (A. 7)  and  (A. 3)  and  expand  term  by  term  in  vector  notation, 

^  1 

—  {T(7)u;}  r(7)u.'  =  -[7ru.'/3  +  7wT  +  f2(-u;)][u;  +  77To»  +  7  X  u;] 

(77  4 

=  —  [( 7  •  u})u>  +  7(0;  •  ui)  —  u  X  u> 

+  (7  •  01)7(7  •  u>)  +  7(7  •  u/)2  -  u>  X  7(7  •  u;) 

+(7  ■  u/)(7  Xw)-f7(u)'(7  x  u;))  —  u>  x  (7  x  u;)]. 

Recalling  the  basic  vector  relations; 

a'-(7xw)  =  7'(u;Xu;)  —  0, 
u ;  x  (7  X  o’)  =  (uj  ■  —  (u>  •  7)0;, 

then  the  right  hand  side  of  (A. 12)  simplifies  to: 

=  j{2(7-«f7  +  2(7-w)[w  +  7Xw]} 

=  k(lT ^)[ilT ^>)l  +  <*’  +  7  x 
=  7Ta;r(7)u;. 


(A. 7) 
(A. 8) 
(4).  To 

(A. 9) 

(A. 10) 

(A. 11) 
(A. 12) 


(A. 13) 

(A. 14) 
(A. 15) 

□ 
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B  Supporting  Computations  for  area  moment  of  inertia  and  B- 
spline  model  for  3-axis  slewing  with  tripod  shaped  appendage 

In  this  appendix  we  summarize  the  computation  of  the  area  moments  of  inertia  of  the  tripod 
metering  truss  model  and  the  computation  by  explicit  integration  of  the  required  matrix 
parameters  obtained  from  the  B-spline  finite  element  model  for  3-axis  slewing  of  a  rigid 
body  with  attached  appendage  consisting  of  a  tripod  with  equilateral  triangular  base.  Only 
those  matrix  parameters  that  are  a  function  of  the  area  moments  are  presented  here.  The 
supporting  computations  for  the  matrix  parameters  that  are  not  a  function  of  the  area 
moments  can  be  found  in  [BBKA88]. 

Computation  of  area  moments  As  described  earlier,  the  metering  truss  was  modeled 
as  a  tripod  of  tubes  located  at  the  vertices  of  an  equilateral  triangle.  A  cross-section  of  the 
structure  appears  in  Figure  2.1.  The  area  moments  of  inertia  can  be  computed  using  the 
following  relations: 


Cr  3 1 ar  +  /x||axij 

(B.l) 

ly  ~  %Icir  4"  -fyllojfjj 

(B.2) 

lz  =  I x  +  /y, 

(B.3) 

where  7cir  is  the  area  moment  of  inertia  of  a  circle  and  is  given  by 

nr4 

Icir  =  “4“’ 

(B.4) 

and  Ix\\axi,  ,  /y||axii  are  additional  terms  that  are  a  consequence  of  the  axis  of  rotation  being 
translated  from  an  axis  through  the  centroid  of  the  circle  to  tht  x  or  y  axis  respectively. 


Figure  2.1:  Cross-section  of  tripod  metering  truss. 
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Beginning  with  the  x-axis,  the  area  moment  associated  with  the  axis  translation  can  be 
computed  as, 


=  nr2d2  +  2 


irr2d2 


From  the  truss  geometry  it  can  be  shown  that, 


d  =  b(  1  -  y  ), 


(B.5) 


( B  .6) 


where  b  is  the  distance  from  the  center  of  the  tripod  base  to  the  center  of  one  of  the  tubes, 
and  /  is  the  distance  from  the  base  of  the  tripod  to  the  end  of  the  truss.  Using  (B.6),  (B.5) 
becomes, 


I xMaxia  (  1  ,  )'  T 


7 rr 


2b2 


r  '  2 

3 

+ 


■d  -  t)2 


3  ,  2  3 irb2  3tt r2b2  2 

=  rr'b-~r  +  ^:- 


(B.7) 

(B.8) 


Finally,  substituting  (B.4)  and  (B.8)  into  (B.l)  gives. 


Ix  = 


37rr4  Zirr2b2  Znr2b2  Zir  r2b2 


+ 


l 


■z  + 


2 12 


(B-9) 


Next,  the  area  moment  of  inertia  for  the  y-axis  is  computed  similiarly  as  the  sum  of 
the  moments  of  inertia  of  three  circles  about  their  centroids  and  the  additional  moments  of 
inertia  of  two  circles  whose  rotational  axes  have  been  translated  (see  Figure  2.1).  The  latter 
term  is  given  by, 

Zirr2d2 

Iy\\axis  —  Z  >  (B.10) 


substituting  (B.6)  into  (B.10)  and  simplifying  we  obtain, 

,  Ztrr2b2  Znr2b2  Znr2b2 


‘y||a 


-2  + 


2  l  2 12 

Finally,  Iy  is  computed  from  (B.2)  using  (B.ll)  and  (B.4)  as, 


(B  11) 


4 


37rr4  3tt  r2b2  37rr262  37rr262 


+ 


/ 


■2  + 


2 12 


(B.12) 


The  area  moment  of  inertia  for  the  z-axis  is  obtained  by  substituting  (B.8)  and  (B.ll) 
into  (B.3)  which  gives, 


I,  = 


37T  r4 


+  Znr2b2  - 


6nr2b2  Znr2b2 


l 


•2  + 


l2 


( B.  13) 
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Computation  of  matrix  FEM  parameters  In  [BBKA88]  the  flexible  appendage  was 
modeled  as  a  beam  with  uniform  cross-section  and  with  corresponding  constant  area  mo¬ 
ments  of  inertia  along  its  length.  However,  in  this  report,  the  tripod  shaped  appendage  has 
area  moments  of  inertia  that  vary  along  its  length  according  to  (B.9),  (B.12),  and  (B.13). 
Therefore,  it  is  necessary  to  recompute  the  four  matrix  parameter  relations  which  are  a 
function  of  the  area  moments  using  the  new  relations.  The  affected  integral  equations  are: 


N„„ 

=  f  {:’pA  +  pHz)}J-, 

Jo 

(B.14) 

Ne 

Jo 

(B.15) 

Nt 

=  [( pI(z)$T(z)  dz, 

Jo 

(B.  16) 

Ke 

=  J  {EI{z)$z(z)$Tz(z)  +  KGdf(2)$r(:)}  dz. 

(BIT) 

Computation  of  N„h 

According  to  (B.9)  and  (B.12)  Ix  =  Iy.  Therefore,  let 

I(z)  =  Ix  =  Iy  =  C\  +  c2~  4-  C3;2  (B.18) 

where 

3tc4  3tt  r2b2 

Cl  "  4  +_ 2~' 

—  3tt  r2b2 

°2  =  }  ’ 

3nr2b2 

C3  “  2 12  ' 

Substituting  (B.18)  into  (B.14)  gives, 

NUh  =  [  {z2pA  +  p(cj  4-  c2r  +  C3Z2)}  dz.  (B.19) 

Jo 


By  direct  integration  of  (B.19)  we  obtain, 

l3  l 2 

=  p[  — (A  4-  c3)  -f  —  c2  4-  lei]. 

Computation  of  Ng 

By  definition  the  right  hand  side  of  (B.15)  has  elements  given  in  terms  of  linear  B-splines 
as 

f(  I(z)B,1_1(z)Bj_1(z)  dz 
Jo 

for  i,j  =  1, . . . ,  N.  Given  the  “hat”  shaped,  linear  B-spline, 

J3‘(z)  :=(?)[=-  ;(*)]B”(.-)  4  (?)[(i  +  2)4 


(B.20) 
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where  the  zero  order  B-spline  is  chosen  with  continuity  from  the  left;  i.e., 

D0/  ,  J  1,  if  Ti  <  x  <  xi+1 
S'(~r):=\0,  else 

Substituting  (B.18)  into  (B.20)  we  obtain, 

—  P  [  (Cl  +  C2Z  +  C3Z2)Bj_1(z)Bj_l{z)  dz 
Jo 

Then  we  have  three  nonzero  cases  to  consider. 

Case  1:  i  =  j  =  N.  Then 

=  pf  (C!  +  c2z  +  c3z2)(^)2{zN  -  ~)2  dz 

J  ’.v-l 

/  l2  l 3 

=  +  C2=N  +  C3~'v2)  ~~  4A^(°2  +  2c3=N’  +  5 Jpc^- 

Case  2:  i  =  j  ±  1.  Here  we  obtain  a  simple  form  for  the  integral  by  transformation  of  the 
variable  and  limits  of  integration; 

|M]y  = 


Case  3:  i  —  j  <  N. 

Mv  = 


Computation  of  NUh$ 

Under  the  assumptions  of  the  3-axis  slewing  model,  the  1  x  N  coupling  coefficient  matrix 
NWhe  has  elements 

=  p  [  (cj  +  c2z  -f  C3r:)B,'(r)  dz , 

Jo 


^  Jo  (Cl  ^  C2C  C3z2}Bi-i(z)B}-i(z)  dz 

P  f  (ci  +  c2z  -f  c3z2)(^  )2(z  -  )( -y  -  z)  dz 

p(t)2  [N[ic  1  +C2(zj  -  x)  +  c3(zj  -  x)2)(£  -  x)x)  dx 
Jo 

l  l2  l3 

^6iV^Cl  +  C2'j  +  C3'^  ~  +  2c3“>)  +  ^T74c3]- 


12jV: 


20  ;V3 


P  [  (ci  +  c2~  +  c3~ 2 ) B}_  1(z)Bj_1(z)  dz 
Jo 

P  f  (C1  +  C2Z  +  C3;2)(  y)2(r  -  Zj_  2)2  dz 

Jzl-  2 

+P  [  {Ci  +  c2z  +  c3z?)(*)2{z}  -  z)2  dz 
J-J-l 

u 

P  /,N(C1  +  c2 (~j  -  *)  +  C3(-J  -  a:)2)(  t)2(2^?  -  i)2dx 

+p  [N(cl  +  c2{zj  -  x)  +  c3(r>  -  x)2)(*)2x2  dx 
Jo 

21  2 12  III3 

+  CiSj  +  C3:^  ~  ZN2^2  +  2c3Cj)  +  1$N*C^ 
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for  i  =  1 . N.  There  are  two  cases  of  interest. 

Case  1:  i  <  N.  Then 

J  (ci  +  C2Z  +  CzZ2)B] (z)dz  —  (^)(ci  +  c2z  -f  c3r2)[z  —  i(£)jdz 

+  r+2m(c!  +  c2z  +  c3z2)[(i  +  2)(£)  -  z]  dz. 

•'“t+t 

Changing  limits  of  integration  by  the  transformation  x  =  z  -  Zi  in  the  first  integral  gives 
J  (  ^)(ci  +  c2~  +  c3c")[;  —  i(ji)}  dz 

=  (t)  fN[c  1  +  c2(x  +  z«)  +  C3(x  +  z,)2]x  dx 
Jo 

l  l 2  / 3 

=  2^V(Cl  +  C2“’'  +  C3r,2)  +  3iV2(C2  +  2c3~*)  +  ^C3 

Changing  variables  (and  limits)  of  integration  in  the  second  integral  yields 

/  (  t  )(ci  +  c2z  -f  c3r2)[(i  +  2)(  £ )  -  r]  dz 

.u 

=  (t)  [”{c i  +  c2(x  +  Zi)  +  c3(x  +  rt)2](2£  -  x)  dx 

l  2 12  III3 

=  2iV(Cl  +  C2Z'  +  C3C,2)  +  3 iV2(C2  +  2c3“l)  +  l2iY3C3' 

From  which  we  can  obtain 

r(  l  p  7/3 

/o  (ci  +  c2=  +  C3'2)SI1(r)  dz  =  —  (ci  +  c2zi  +  c3r2)  +  — (c 2  +  2c3Zi)  +  ^jc3. 

Case  2:  i  =  A". 

[  {c\  +  c2z  +  c3z2)Bi(z)dz  =  f  (cj  +  c2r  +  c3r2)( ^ )[r  —  r^_x]dr 

*0  J:N-i 

=  j  t  (t)[ci  +  c2(x  +  zN)  +  c3(x  +  z^)2](x  +  £)  dx 

J~Tr 

l  l2  l3 

=  -^(Cl  +  C2ZiV  +  c3z^)  +  — ^(c2  +  2  C3Z,v)  +  ^JC3 

Computation  of  Kg 

Under  the  assumptions  of  the  3-axis  slewing  model,  the  N  x  N  coefficient  matrix  Kg  has 
elements 

\K,]ii  =  j'{EI(z)<i>:(z)*Tz(z)  +  KGA*(z)*T(z)}dz, 

=  EIN"  +  kGAN. 
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where  N"  has  elements 


=  l‘  £s,'.1(.-)fl‘.,(  =  )(c1  +  c,z  +  c,z2)  d = 


and  N  is  given  in  [BBKA88].  For  the  N  x  N  tridiagonal  matrix  N"  and  i,j  =  1, . . . ,  N  we 
obtain  three  nonzero  cases. 

Case  1:  i  =  j  <  N . 


Case  2:  i  =  j  =  N. 


f  (^)2(ci  +  c2 z  -f  C3Z2)  dz  +  f  (^)2(cj  +  C2 z  +  C3Z2)  dz 
Jo  Jo 

2  N  21 

~TC  1  +  c2  +  -ZT7C3. 


[N'%  =  [N(t)2(ci  +  C2r  +  C3C2)  dz 

Jo 

N  1  l 

-  JCl  +  2C,  +  3NC^ 


Case  3:  i  =  j  ±  1.  For  i  —  j  +  1  we  obtain. 


[N  ]«'j  —  [  ~ ( ^)‘(ci  +  c2~  +  c3z2)  dz 

J  -1  +  1 

=  [,  ~{t)2(ci  +  C2(x  +  Zi)  +  c3(x  +  zt)2)  dz 

J  N 

=  -y(°l  +  C2=i  +  C3~2)  -  -(c2  +  2c3Ci)  -  ^C3. 


Similarly,  for  i  =  j  —  1 


TSI-89-12-12-WB 
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